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ABSTRACT 

Context. Type la Supernovae (SNe la) have been widely used in cosmology as distance indicators. However, to fully exploit their 
potential in cosmology, a better control over systematic uncertainties is required. Some of the uncertainties are related to the unknown 
nature of the SN la progenitors. 

Aims. We aim to test the use of integral field unit (IFU) spectroscopy for correlating the properties of nearby SNe la with the properties 
of their host galaxies at the location of the SNe. The results are to be compared with those obtained from an analysis of the total host 
spectrum. The goal is to explore this path of constraining the nature of the SN la progenitors and further improve the use of SNe la in 
cosmology. 

Methods. We used the wide-field IFU spectrograph PMAS/PPAK at the 3.5m telescope of Calar Alto Observatory to observe six 
nearby spiral galaxies that hosted SNe la. Spatially resolved 2D maps of the properties of the ionized gas and the stellar populations 
were derived. 

Results. Five of the observed galaxies have an ongoing star formation rate of 1-5 M Q yr _1 and mean stellar population ages ~ 5 Gyr. 
The sixth galaxy shows no star formation and has an about 12 Gyr old stellar population. All galaxies have stellar masses larger than 
2 x 10 10 M o and metallicities above solar. Four galaxies show negative radial metallicity gradients of the ionized gas up to -0.058 
dexkpc -1 and one has nearly uniform metallicity with a possible shallow positive slope. The stellar components show shallower 
negative metallicity gradients up to -0.03 dexkpc -1 . We find no clear correlation between the properties of the galaxy and those of 
the supernovae, which may be because of the small ranges spanned by the galaxy parameters. However, we note that the Hubble 
residuals are on average positive while negative Hubble residuals are expected for SNe la in massive hosts such as the galaxies in our 
sample. 

Conclusions. The IFU spectroscopy on 4-m telescopes is a viable technique for studying host galaxies of nearby SNe la. It allows one 
to correlate the supernova properties with the properties of their host galaxies at the projected positions of the supernovae. Our current 
sample of six galaxies is too small to draw conclusions about the SN la progenitors or correlations with the galaxy properties, but the 
ongoing CALIFA IFU survey will provide a solid basis to exploit this technique more and improve our understanding of SNe la as 
cosmological standard candles. 

Key words. Galaxies: general - Galaxies: abundances - (Stars:) supernovae: general 



1. Introduction 

More than a decade ago the observations of Type la Supernovae 
(SNe la) led to the discovery of the accelerating expansion of the 
Universe and the need for an unknown repulsive force to drive it 
( [Perlmutter et al.|1999||Riess et al.|1998| ). Understanding the na- 
ture of this force - now dubbed "dark energy" - is an outstanding 
goal of astrophysics and cosmology. It is now well-understood 
that no sing le observational techn ique will be able to achieve this 
alone (e.g. Albrecht et al. 2006). The combined constraints of 
many techniques will be needed to measure the equation-of- state 
parameter of the dark energy (the ratio between the pressure and 
the density) and eventually its evolution over the cosmic time. At 
present, SNe la is the most matured and well-understood tech- 
nique to accurately trace the cosmic expansion history and will 
contin ue to play an essential role in future cosmological experi- 
ments ( [Albrecht et al.|200"6 ). However, the SNe la technique is 
affected by several systematic uncertainties, which need to be 



* Based on observations collected at the Centro Astronomico 
Hispano Aleman (CAHA) at Calar Alto, operated jointly by the Max- 
Planck Institut fur Astronomie and the Instituto de Astroffsica de 
Andalucfa (CSIC) 



reduced to a level below ~ 2% to differentiate between different 
dark energy models. 

The use of SNe la in Cosmology relies on the empirically 
established tight relation between their light curve width and 
peak luminosity, which allows one to measure the luminosity 
distance with an accuracy of ~ 7% (Phillip set al.|1999] ), and on 
the assumption that the (standardized) peak luminosity of SNe 
la does not change over the cosmic time. There is now obser- 
vational evidence that the slope of the "light curve shape - peak 
luminosity" relati on does not depend on the redshift or the host 
galaxy mass ([Conley et al.|2 01 1 ; Sulliva net al.|201l| ). However, 
Sull ivan et al. 1(2010) and |Kelly et al.|p010| ) havefound that the 
offsets of the SN la peak magnitudes from the best- fitting Hubble 
line (from now on Hubble residuals or HR) correlate with the 
host stellar mass. Together with the mass-metallicity relation for 
galaxies (e.g., Tremon tiet al.|2004] ) and the overall increase of 
the metal content of the universe with the cosmic time, this may 
be an indication for possible luminosity evolution of SNe la at a 
level of 0.05-0. 10 mag. 

It is now generally agreed that SNe la are the result 
of thermonuclear disruption of carbon/oxygen (C/O) white 
dwarfs (WD), which ignite explosively when they approach the 
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Chandrasekhar limit M ch ~ 1.38M© ( [Hoyle & Fowler |[196Q| ). 
However, there has been little observational evidence of the 
exact evolutionary scenario that leads to the explosion. C/O 
WDs are the end product of the evolution of stars with masses 



1.5 - 7M Q (e.g., see |Becker & Iben|1980[ 

ZDs is - 



[T999]). The upper mass limit for C/O W 



Dominguez et al. 
1.1 Mq (e.g., see 



DommguezeFaT][T999| |Salaris et al.pX)9| |Weidemann|[T987l ) 



and therefore a mechanism that allows the WD to gain addi- 
tional mass of at least ~ 0.3 - OAMq is needed. In the single 
degenerate (SD) scenario the WD accretes m ass from a non- 
degenerate companion star in a binary system ( Whela n & Iben| 
|1973| ). However, the exact physical mechanism of the WD mass 
growth has not yet been identified. In the double-degenerate 
(DD) scenario two C/O WDs in a binary merge after loosing or- 
bital angular momentum by gravitational wave ra diation pben &| 
|Tutukov| 1984j |Paczynskip985l |^bbink|1984| ). While consid- 
ered the most viable, both scenarios have considerable uncertain- 
ties (see, e.g. Hillebrandt & Niemeyer 2000; Maoz & Man nucci] 
|20TTT ). 

The numerical simulations of thermonuclear SN la explo- 
sions suggest that the properties of the exploding WD may sig- 
nificantly influence the peak luminosity, the "light curve width - 
lumi nosity" relation a nd colors of the resulting sup ernovae (e.g., 
Bravo et al.|[20T0| JDormnguez et al.|[200T| |Hoflich et al.p 998; 



Kasen et al.l|2009[ |Ropke et al.| |2006; Um eda et al.||1999| ). On 



the other hand, the properties of the WD just before the ignition 
(the central density, metallicity and C/O ratio) are sensitive to the 
properties and the evolution of its progenitor binary star, and to 
the subsequent WD mass growth mechanism. For example, the 
SD channel can produce Mch WDs with slightly different struc- 
ture and chemical composition depending on the mass of the 
WD at the moment when the accretion started (e.g., Dominguez 
et al.|200T] ). In the DD scenario, the outcome of the merger may 



depend on the mass ratio of the two WDs. In addition, one may 
expect that some properties of SN la progenitor stars will evolve 
with cosmic time, e.g. metallicity. Therefore, possible evolution 
of the properties of SN la progenitors or, if more than one evo- 
lutionary channel exist, evolution of their relative contribution 
to the SNe la population, may introduce systematic uncertain- 
ties in SN la cosmology and potentially bias the cosmological 
results from the future large SN surveys (e.g., see Nordin et al. 
|2008l|Sarkaretal.|2008| ). 

To date no progenitor of a SN la has been unambiguously 
identified and/or observed and information about the SNe la pro- 
genitors has been inferred indir ectly. Mannuc ci et al. (2005 ) and 
Scannapieco & Bildsten (2005) studied the SN la rate as a func- 



tion of redshift and host galaxies properties. Both studies found 
that the SN la rate depends on both the on-going star forma- 
tion rate (SFR) and total galaxy stellar mass. This result was 
also confirmed by other s ([Dahlen et al.||2008[ |Mannucci et al. 



2006 l|Maozetal.||20TTl |Neill et al.||2006[ |Pritchet eTaLl |2008 ; 



Sulli van et al.|2006) and appears to suggest that at least part of 
SNe la are associated with the young stellar population capable 
of producing SN la with shor t delay time < 400 Myr. |Totani| 
|et al.| ( |2008| ) and |Maoz et al.] ( [20101 ) have shown that the de- 
lay times from star formation to SN la explosions between the 
shortest time probed < 400 Myr and 10 Gyr are distributed as a 

power law with slope 1 . This delay time distribution (DTD) 

strongly favors the DD scenario. The SD scenario may also ex- 
plain this DTD (Hac hisu et al.||2008] ) but the efficiency of the 
symbiotic channe l (WD+re d giant) needs t o be significantly in- 
creased (see, e.g., Maoz & Mannucci 201 1 ). The early discovery 
of SN 201 lfe in M101, the nearest SN la in 25 years, provided 
the first real possibility to constrain the properties of a progen- 



itor star of an SN la ( [Bloom et all2012| |Chomiuk et aLl|20T2l 
Nugent et al. 201 1]). The results reinforce the conclusion that the 
exploding star is a C/O WD and seem to rule out all but a de- 
generate star as its companion, thus favoring the DD scenario. 
On the other hand, based on high-resolution spectroscopy of a 
sample of nearby SNe la [Sternberg et al.| ( |2011| ) favor the SD 
scenario. 

Many studies have shown that the intrinsically luminous SNe 
tend to occur in star-formin g hosts, while the faint SNe pre- 
fer passive ones (e.g., [Brandt et al.|2010||Gallagher et al.|2005 



2008|| Hamuy et al.|19 96, 2000; |Howell et _ al.|2009[|KelIy et al 



2010; Neill et al. 2009; Raskin et al. 2009 HSullivan et al.|2 010, 
2006). Galla gher et al.| p008| ) found that the Hubble residuals 



correlate with the global host metallicity. |Sullivan et al. 



and |Kelly et al. ( 2010|) found such a corr elation with t 
stellar mass. However, 



(2010) 



le host 

iowell et aL]([2009 ), who used the galaxy 



stellar mass as a proxy for the metallicity, found no such correla- 
tion and suggested that instead the progenitor age may be a more 
important parameter. 

All studies of SN la hosts galaxies conducted so far, except 
that by Rask in et al.| ( |2009| ), were based on an analysis of the 
global photometric or spectroscopic properties of the host galax- 
ies. In this paper we take a different approach. We use for the 
first time integral field unit (IFU) spectroscopy at intermediate 
spectral resolution to study a sample of host galaxies of local 
SNe la (z ~ 0.02). This approach has an advantage over the pre- 
vious studies because it allows us to derive spatially-resolved 
two-dimensional (2D) maps of host galaxy properties, e.g. the 
heavy element abundance in the interstellar medium (ISM). The 
intermediate spectral resolution makes it also possible to use 
full-spectrum fitting techniques to derive 2D maps of the prop- 
erties of the stellar populations. The main objective of this pilot 
work is to test the methodology to correlate the properties of the 
SNe with the properties of the gas and the stellar populations at 
the location of the SN explosion, in addition to the global host 
properties. By analyzing the properties of the stellar populations 
we also aim to constrain the nature of the SNe la progenitors. 

Throughout the paper we assume the concordance cosmo- 
logical model with £l M = 0.27, Da = 0.73, w = -1 and 
= 0.71. 



2. Observations and data reduction 

2.1. Target selection 

The list of targets for this program was selected from a sample 
of spiral galaxies that hosted SNe la for which the important pa- 
rameters such as luminosity, extinction, intrinsic color indices, 
luminosity decline rate, and deviation from the Hubble diagram 
have been accurately measured. The galaxies were carefully ex- 
amined and selected to fulfill four additional criteria: 

1. to have angular size ^40-60 arcsec; 

2. to be nearly face-on; 

3. the SN lies on a high surface brightness location in the 
galaxy; 

4. be observable at airmass less that 1.3 to minimize the effect 
of the differential atmospheric refraction. 

The first two requirements maximize the use of the large 
field-of-view (FOV) of the IFU instrument and minimize the pro- 
jection effects when correlating the SN and its local host galaxy 
properties. The third requirement ensures that we will obtain a 
good signal-to-noise ratio (S/N) of the spectra at the location 
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Table 1. Supernovae and details of their host galaxies: morphological type, Milky Way dust reddening, offsets from the host nucleus, 
de-projected galactocentric distance, inclination, and position angle. 



SN 


Host 


Type a 


z h 


E(B - V)mw c 


RA offset 
[arcsec] 


DEC offset 
[arcsec] 


DGD^ 
[kpc] 


i b 
[deg] 


PA* 7 
[deg] 


1999dq 


NGC 976 


SAb 


0.0144 


0.110 


-4.0 


-6.0 


2.4 


36.5 


77.6 


1999ej 


NGC 495 


SB(s)0/a 


0.0137 


0.072 


+ 18.0 


-20.0 


7.7 


44.0 


47.6 


2001fe 


UGC5129 


SBbc 


0.0134 


0.022 


-13.5 


-0.1 


3.9 


46.0 


14.1 


2006te 


CGCG 207-042 


SB(r)bc 


0.0315 


0.046 


-5.5 


-1.7 


4.9 


39.4 


69.9 


2007A 


NGC 105 NED02 


SBc 


0.0173 


0.073 


-1.2 


+ 10.1 


3.6 


36.4 


77.3 


1997cw 


NGC 105 NED02 


SBc 


0.0173 


0.073 


+8.0 


+4.0 


4.0 


36.4 


77.3 


2007R 


UGC 4008 NED01 


SAa 


0.0308 


0.047 


-1.9 


-3.9 


3.4 


48.1 


76.4 



(a) based on SDSS pseudo-color images. ^ this work. Derived from analysis of the Ha velocity field, except for the host of SN 1999ej, for which 
the stellar velocity field was used. (c) from Schlegel et al. ( 1998 1. 



of the SN, and in particular will allow us to access the absorp- 
tion and emission line spectra. From the large sample of galaxies 
six objects were observed. The galaxy details and the SNe off- 
set from the nucleus are given in Table[T] All SNe are normal 
SNe la, except for SN 1999dq, which was classified as a pecu- 
liar 1991T-like event ( |Jha et al.|1999] ). Tableggives the SALT2 
x\ and C parameters of the supernova light curves ( |Guy et aL] 
12007) ) taken from | Amanullah et al.| ( |2010| ), the offset from the 
best-fit Hubble line A//, and the AM15 parameter, which shows 
how much the SN 5-band magnitude has declined during the 
first 15 days after the time of the B band maximum. x\ and C 
are parameters related to the SN light curve shape and B - V 
color index at maximum, respectively. They are used to stan- 
dardize the observed 5-band peak magnitude .Sobs via the rela- 
tion #std = #obs + ax\ - J3C , with a - 0.121 and J3 = 2.51 as 
per I Amanullah et aL]p010| ). A// was computed after first cor- 
recting the redshifts of the SNe for large-scale coherent galaxy 
motions in the local universe based on the models of Hudson 



et al. (2004). The accuracy of this correction is estimated to be 
~ 150 km s" 1 and a random peculiar velocity of 150 km s" 1 is 
added to the uncertainty of A//. 

It should be noted that for this pilot project the selection cri- 
teria are solely optimized to maximize the quality of the obser- 
vations and facilitate the analysis. We focus on late-type galaxies 
because one of our goals is to correlate the SN properties with 
the properties of the ISM determined from the ionized gas. This 
leads to strong biases, however, e.g. the galaxies in our sample 
are bright, massive, and likely metal-rich. 

2.2. Observations 

The six galaxies were observed on November 14 and 15, 2009 
at the 3.5m telescope of the Calar Alto observatory using the 
Potsdam Multi- Apert ure Spectrograph (PMAS,|Roth et al |2005| ) 
in the PPAK mode JKelz et al.||2506| |Verheijen et al.||20Q4| ). 
The atmospheric conditions were variable with occasional thin 
clouds interrupting the observations. The seeing varied between 
1.5"and2.2". 

The PMAS instrument is equipped with a 4kx4k E2V#231 
CCD. We used a set-up with the 600 lines mm -1 grating V600 
and 2x2 binned CCD, which provided a wavelength range of 
-3700-7000A with a spectral resolution of ~ 5.5A. The PPAK 
fiber bundle of PMAS consists of 382 fibers with 2.7"diameter 
each, 331 of which (science fibers) are ordered in a single hexag- 
onal bundle that covers a FOV of 72 // x64 // . Thirty-six additional 
fibers form six mini-bundles (sky -bundles), which are evenly 
distributed along a circle of ~ 9 0"radius and face the edges of 
the central hexagon (see Fig. 5 in Kelz et al.|2006 ). The remain- 



ing 15 fibers are used for calibration and can only be illuminated 
with the PMAS internal calibration unit. For a detailed descrip- 
tion of the PPAK fiber bundle we refer the reader to IKelz et aLl 
( 2006 ). Some details that are relevant for the data reduction are 
also given in Appendix [A[ 

For each object three 1800-sec long exposures were ob- 
tained. Because the filling factor of a single PPAK exposure is 
-65%, we adopted a dithering pattern with the second and the 
third exposures offset by A(R.A.,Decl.)=(1.56, 0.78) and (1.56, 
-0.78) arcsec with respect to the first exposure to ensure that ev- 
ery point within the FOV was spectroscopically sampled. Before 
and after the science exposures, spectra of HgNe and continuum 
halogen lamp were obtained to wavelength-calibrate and trace 
the spectra. The spectrophotometric standard stars Feige 34 and 
BD+25 3941 were observed to measure the sensitivity function 
of the instrument. In addition, series of exposures of blank sky 
regions were obtained during twilight and were used to equalize 
the fiber- to-fiber throughput variations. 

2.3. Data reduction 

The pre-reduction of the CCD images was performed with 
IRAljjand the rest of the reduction with our own programs writ- 
ten in IDL. Each individual science pointing was reduced inde- 
pendently. After the standard CCD reduction steps of bias sub- 
traction, flat-field correction and removal of cosmic ray hits, the 
spectra were traced, extracted, wavelength- and flux-calibrated, 
and finally sky- subtracted. At the final step the three pointings 
were combined into a final 3D data-cube, taking into account 
the differential atmospheric diffraction. The full details of the 
data reduction are given in Appendix [B] 

Three of the galaxies in our sample also have SDSS spectra. 
This allowed us to check the relative flux calibration of our spec- 
troscopy. Spectra within an aperture of 3" diameter centered on 
the galaxy nucleus were extracted from the data-cubes to emu- 
late the SDSS spectra. The comparison, after our spectra were 
scaled to match the flux level of the SDSS spectra, is shown in 
Fig. [I] It demonstrates that the relative flux calibration of our 
spectra is excellent and matches SDSS to within a few percent. 

The absolute flux scale of the data-cubes was set using the 
SDSS imaging. SDSS g and r magnitudes of the galaxies were 
computed within an aperture of 20" diameter. Spectra within the 
same aperture size were extracted from the data-cubes and syn- 
thetic g and r magnitudes were computed. The g and r scale 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Table 2. SALT2 |Guy et al |2007|) x x and C parameters of the 
SNe from |Amanullah et al.| pOTO| ), the offset from the best-fit 
Hubble line A//, and the AM15 parameter. The uncertainties of 
the parameters are given in the parentheses. 



SN 


Xi 


C 


A/i 


AM 15 < 


1999dq 


0.89 (0.12) 


0.13 (0.01) 


-0.33 (0.09) 


0.96 


1999ej 


-2.08 (0.43) 


0.07 (0.05) 


0.46 (0.20) 


1.49 


2001fe 


0.41 (0.18) 


0.03 (0.02) 


-0.06 (0.09) 


1.03 


2006te 


-0.36 (0.18) 


-0.04 (0.02) 


0.16(0.08) 


1.15 


1997cv/ 


0.79 (0.25) 


0.40 (0.03) 


-0.08 (0.13) 


0.97 


2007A 


-0.04 (0.14) 


0.18(0.01) 


0.15 (0.07) 


1.10 


2007R 


-1.76 (0.16) 


-0.07 (0.02) 


0.23 (0.07) 


1.42 



(a) calculated from x\ with the relation given in |Guy et al. ( 2007| ); the 
first photometric observation was taken ~15 days past maximum and 
the photometric parameters are rather uncertain. This SN was included 
in the analysis because it is in the same host as SN 2007A. 

factors that provided the match of the synthetic magnitudes to 
the observed ones were computed and the average of the two 
was applied to the data-cubes. We note that the g and r scale 
factors coincided to within 3%, which additionally supports our 
conclusion that the relative flux calibration is accurate. 



3. Data analysis 

The individual spectra in the data-cubes were analyzed to de- 
rive 2D maps of the properties of the galaxies. This included 
the properties of the ionized gas and the stellar populations. The 
properties of the galaxies at the SN position and the galaxy cen- 
ter were derived by interpolating the 2D maps. The galaxy cen- 
ters were computed from the data-cubes and the SN positions 
were computed with respect to it, using the offsets quoted in the 

pOQg) . 



Jha et al. 



discovery IAU circulars and 

As previously mentioned, one of the main goals of this study 
is to test the feasibility of using IFU spectroscopy to compare the 
properties of the host as derived from integrated spectroscopy 
to those derived from spatially resolved spectroscopy. For this 
purpose, we also analyzed for each galaxy the total spectrum 
formed by simply summing all spaxel spectra in the data-cube. 
This simulates an observation of the same galaxy with long-slit 
spectroscopy as is performed for high-redshift galaxies. 

For each galaxy the analysis was also performed on az- 
imutally averaged spectra at several de-projected galactocentric 
radii. This was performed as an alternative way to derive the ra- 
dial dependence of the galaxy properties such as the metallicity. 
To compute the azimutally averaged spectra we first computed 
the de-projected galactocentric distance of each spaxel with the 
position angle and inclination (Table [1} compu ted from the anal- 
ysis of the Ha velocity map^](see Sec. |3. 1.2 }. The spectra were 
corrected to rest-frame wavelength with the stellar velocities es- 
timated from the fits to the absorption line spectrum (Sec. 3.2.4 ). 
Finally, for each galaxy the spectra within several (4 to 6) radial 
bins were averaged. Because the stellar velocities were used to 
correct the spectra to rest frame, we used only the spectra that 
had a sufficiently high S/N to allow fitting with STARLIGHT. 

All quoted uncertainties of the derived quantities are statis- 
tical and do not include systematic and intrinsic uncertainties of 
the methods, which will be additionally discussed when appro- 
priate. In the next sections we present the main steps in the data 
analysis. 
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Fig. 1. Comparison between SDSS spectra (red) and spectra ex- 
tracted from our data-cubes within an aperture of 3 " diameter 
centered on the galaxy nuclei (black). 



3.1. Properties of the ionized gas 

The presence of nebular emission lines in the galaxy spectra al- 
lows us to study the properties of the ionized gas such as its oxy- 
gen abundance and ionization state, and to derive other important 
properties such as the star formation rate, dust extinction, etc. 
Some of the methods used to derive these quantities, for example 
the strong line methods to estimate the gas metallicity, can only 
be applied if the ionization source is exclusively arising from the 
stellar radiation. For this reason and to search for possible AGN 
contamination, we used the diagnostic diagram [O m] A5OO1/H0 
vs. [N11] ^6584/Ha (BPT diagram; [BaTdwin et al.|198l) (Fig.g). 
The central spaxels that fall in the AGN area of the diagnostic 
diagram according to the Kauffm ann et al.| ( [2003| ) criterion were 
excluded from the relevant parts of the analysis. 



3.1 .1 . Emission line fluxes 

Five of the six galaxies in our sample show strong neb- 
ular emission lines. The fluxes of the prominent emis- 
sion lines [On] ,13727, H# [Om] /U4959/5007, Ha, and 
[N11] /L16549/6584 were used in the analysis. Whenever possi- 
ble, the [S 11] AA61 16/6731 lines were also measured. In the spec- 
tra of galaxies the emission lines are superimposed on the un- 
derlying stellar absorption spectrum. The stellar absorption lines 
can bias the measurement of the emission line fluxes, an effect 



2 For NGC 495 the stellar velocity map was used because this galaxy 
shows no emission lines. 



that is especially prominent in the Hj3 line (Fig. B.3 ). Therefore, 
to measure the emission line fluxes accurately, the stellar absorp- 
tion spectrum needs to be subtracted first. For this w e used the 
STARLIGHT software ( |Cid Fernandes et al ||2Q05] ). All spec- 
tra that had S/N greater than 5 at ~ 4600A were fitted with 
STARLIGHT and the emission line fluxes were measured on 
the continuum-subtracted spectrum. For the remaining the spec- 
tra the measurements were made without continuum subtraction. 
Each emission line was fitted with a single Gaussian plus a lin- 
ear term, and the area under the Gaussian was taken as flux es- 
timate. Details of the adopted procedure and the Monte Carlo 
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simulations that were performed to estimate the uncertainties of 
the line fluxes are given in Appendix [C] 

3.1.2. Ha velocity field 

The fitted positions of the strongest of all emission lines, Ha, 
provide the best estimate of the gas velocity field. These fields, 
shown in Figs. [3][5] were analyzed with the methods and IDL 
programs developed by Krajnovic et al. ( 2006| ). The program 
analyzes the velocity field at several radii and for each of them 
returns the inclination and position angle and quantifies the de- 
gree of deviation from pure disk rotation. From this analysis we 
also derived the redshift, the average position angle (PA) and in- 
clination i for each galaxy, which are listed in Table [T] 



3.1 .3. Extinction, Ha flux, and star formation rate maps 

For the purpose of the following analysis the measured line 
fluxes were corrected for dust extinction using the observed 
B aimer decrement I(Ha)/I(Hj3) and assuming a foreground dust 
screen. For the intrinsic B aimer decrement I(Ha) I KH/3\ m a 
value of 2.86 was assumed, which is appropriate for case-B re- 
combination with electron t emperature T e - 100 00 K and elec- 
tron density 10 2 c m 3 (e.g.,|Osterbrock & Ferland 2006| . The dust 
is described by the Fitzpatrick (1999 ) law with R v =3.1. 

The extinction-corrected Ha flux was converted into instan- 



taneous SFR using the|Kennicutt|( 1998 ) relation: 



SFR [Af yr" 1 ] = 7.9 x 10" 42 L(Ha), (1) 
where L(Ha) is the Ha luminosity in units of erg s" 1 . 

3.1.4. ISM oxygen abundance 

The most accurate method to measure the ISM abundances 
- the so-called direct method - involves determining of the 
ionized gas electron temperature, T e , which is usually esti- 
mated from the flux ratios of auroral to nebular emission 



lines e.g. [Q in] >U49 59/5007/[O m] /14363 (e.g. jlzotov et al. 
2006; Stasinska 2006). However, the temperature- sensitive lines 
such as [Om]/14363 are very weak and difficult to mea- 
sure, especially in metal-rich environments. A careful exam- 
ination of our data-cubes revealed that the [Om]/l4363 line 
was not present. For this reason we used other strong emis- 
sion line methods to determine the gas oxygen abundance. 
Many such methods have been developed throughout the years, 
the most commonly used being R 2 3 =([On] /13727+[Oiii] 
/U4959/5007)/Hff ratio-based methods dKewley & Dopita|2 002 ; 
Kobulnicky & Kewley|2004| |McGau gh 1991; |Pagel etal. 1979 
Pilyugin||2001| |Pilyugin & Thuan||2005[ |Tremonti et al.pOCH 



2004[ |Storchi-Bergmann et al. 



^5007)/HyS)/([Nn]^6584/Ha)] (Alloin et al. 19791 |Pettini & 



Zarit sky et al.|1994|), N2=log[[Nn] /[65 84/Ha] (Pettini & Pagel 



1994]) and Q3N2=log[([Om] 



Pagel 2004). More recently, |Liang et al.| ( |20Q7[ |2Q06| ) and|Yin 



et al.| (|2007) have verified and re-calibrated these and other 
strong-line methods using Sloan Digital Sky survey (SDSS) 
spectroscopy. 

Unfortunately, there are large systematic differences between 
the methods, which translate into a considerable uncertainty 
in the absolute metallici ty scale (for a recent review see, e.g., 
Kewley & Ellison 2008 ). In particular, there is ~ 0.4 dex dif- 
ference between the so-called empirical and theoretical strong- 
line methods. The empirical methods are calibrated against H II 
regions and galaxies whose metallicities have been previously 



determined b y the direct method, e.g. Q3N2 and N 2 (Pettini 
& Pagel|[2004l ), R 23 - P ( jPilyugin & Thuan|[2005] ). The the- 
oretical methods, on the other hand, are calibrated by match- 
ing the observed line fluxes with those predicted by theoretical 
photoionization models (most of the R23 -based methods, e.g., 
Kewley & Dopita|2002[ |Kobulnicky & Kewley |2004[|McGaugh| 



1991[ |Tremonti et al.||2004 ). The cau se of these discrepancies 
is still not well-understood. Recently [Moustakas et al.| d2010 ) 
discussed this problem and concluded that the empirical meth- 
ods may underestimate the m etallicity by a few tenths of dex 
(see also [Peimbert et al. 2007 ), while the theoretical methods 
overestimate it. In this situation, we followed the recommenda- 
tion of Kewley & Ellison (2008| to use one method to compute 
the metallicities in all galaxies and discuss the results in relative 
sense, and use another method to confirm the observed trends. 
As our primary metho d we used the empirical 03N2 method of 
Pettini & Pagel ( 2004 ) (PP04 from now on) and checked the re- 



sults with the theoretical R23 method of Kobulnicky & Kewley 



(2004) (KK04 from now on). Both methods have advantages 
and disadvantages, which have been discussed in several papers 
(e.g., |Kewley & ElUson|2008||Yin et al.|2007) . 



3.2. Stellar populations 

The star formation history and chemical evolution of a galaxy 
is imprinted in the properties of its present-day stellar popula- 
tions. Determining the properties of the stellar populations in 
the galaxies has been a major research topic in astrophysics and 
through the years many different methods have been used, rang- 
ing from analysis of the color-magnitude diagrams (CMD, [Faber| 
1972) to equivalent widths of absorption lines (e.g., the Lick 
indices, [Worthey et al.||1994] ). However, in most galaxies sev- 
eral stellar population are simultaneously present. Disentangling 
their contribution to the galaxy spectrum is a very difficult task 
because of various astrophysical and numerical degeneracies. 



3.2.1 . Full-spectrum fitting technique 

Recently, the so-called evolutionary population synthesis meth- 
ods dBruzual & Charlot||2003l |Tinsley|[T968l |Vazdekis et aL 
1996b coupled with full-spectrum fitting techniques (e.g, |Cid 



Fernande set al.j2005l|Koleva et al.|2 008 ; MacArth uretaLl2 009 ; 
|Reichardt et al.|200T] [Vazdekis & Arimo tofl999| ) have emerged 
as powerful means to analyze galaxy spectra. The evolutionary 
population synthesis methods produce synthetic galaxy spectra 
using as input theoretical evolutionary tracks, libraries of stellar 
spectra, initial mass function (IMF), and prescriptions for star 
formation and chemical evolution. The models are then com- 
pared to the observed spectra to infer the properties of the stel- 
lar populations that contribute to the formation of the observed 
spectrum. One possible approach is to fit the observed spec- 
trum with a linear combination of model spectra of single stel- 
lar populations (SSP) of different ages and metallicities (e.g., 
Cid Fern andes et al.|2005[ |Koleva et al.|2008[ |MacArthur etaEl 
2009 ). The fitting returns the contribution of the different SSPs 
(called population vector) that best describe the observed spec- 
trum, which then can be used to study the stellar populations 
of the galaxy. However, because of astrophysical and numer- 
ical degeneracies, and the presence of noise in the observed 
spectra, it is well-known that the solution may not be unique 
and the results should be interpreted with caution (e.g, see the 
discussion in |Cid F ernandes et al. 2005 ). The best known is 
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the age-metallicity degenerac)^] where young metal-rich stel- 
lar populations are confused with older metal-poor ones (see 
for example Fig. 10 in|Bruzual & C hariot 2003 ). As noted by 
Bruzual & Chariot ( |2003 ), while the shape of the stellar con- 
tinuum is roughly the same, the strength of the metal lines in- 
creases. Therefore, analyzing well-calibrated spectra with high 
S/N and spectral resolution to resolve the absorption lines has 
the potential to brake the age-metallicity degeneracy. In addition, 
uncertainties in the input ingredients needed for computing the 
SSPs, such as non-uniform coverage of the age/metallicity pa- 
rameter space of the stellar libraries, IMF and the difficulties in 
describing some phases of the stellar evolution (e.g., the thermal- 
pulsating asymptotic giant branch phases), add even more un- 
certainties when interpreting the results (see, e.g., |Conroy et al. 
^009] ). 



3.2.2. Choice of the base 

In this study we used the STARLIGHT code described in |Cid| 
Fernande s et al.| ( |2005| ) and |Asari et aL] ( |2007| ) coupled with a 
version of the Bruzual & Chariot SSP models based on the new 
MILES spectral library ( Sanchez-Blazquez et al.|2006 ). The se- 
lection of the SSP basis is important for any full- spectrum fitting 
algorithm and the interpretation of the results. To minimize the 
computing time one should select few SSPs that are maximally 
independent and at the same time are capable of reproducing the 
variability of the full SSP set for a given metallicity. If a large 
basis is selected, many of its components will be close neigh- 
bors. This will lead to increased non-uniqueness of the solution 
and increase the time for the fitting algorithm to converge. On 
the other hand, if too small a basis is selected, it will not be able 
to capture the full variability of the SSP models, the fits may be 
poor, and the results will be unreliable. 

In our work we used the following approach to select the ba- 
sis. For a given metallicity all SSPs were normalized to the flux 
in the 4600-4800 A interval. Then the evolution of the flux in 
seven spectral windows in the range 3700-7000 A was tracked as 
a function of the SSP age. The goal was to identify age intervals 
where the flux in all seven spectral windows evolves linearly (or 
close to) with time. If such intervals exist, then the SSPs within 
them are not independent; all SSPs in a given interval can be 
closely reproduced as a linear combination of the two SSPs at 
the extremes. By selecting the basis at the ages connecting the 
linear intervals we form a small independent set of basis vectors, 
which at the same time can reproduce the SSPs at all other ages. 
Following this approach we were able to select Af* = 16 or 17 
SSPs per metallicity that formed our fitting basis of 66 SSP mod- 
els with ages between 1 Myr and 1 8 Gyr, and four metallicities 
Z=0.004, 0.008, 0.02 (the solar metallicity) and 0.05. 



3.2.3. Voronoi binning 

To increase the S/N in the outer parts of the galaxies the data 
cubes were spatially binned using adaptive V oronoi tessellations 
( [Cappellari & Copin|2003l |Diehl & Statler|2006| ). The binning 
of the spaxels was determined from the S/N measured in the 
interval 4580-4640 A, after discarding the spectra with S/N<1. 
The targeted S/N of the binned spaxels was S/N-20, with the 
exception of the host of SN 2006te, for which a lower S/N of 



3 Dust reddening also adds to this problem, partly because the dust 
extinction law may be different in different galaxies; galaxies in the 
Local Group are a good example for this. 

4 circa 2007; unpublished 



15 was used. To keep the spatial resolution reasonably small, 
an upper limit of the size of the bins was also imposed: 5 for 
NGC976, 17 for NGC495, and 12 for the remaining four. 

3.2.4. STARLIGHT fits 

The Voronoi-binned spectra along with the un-binned ones were 
fitted with the STARLIGHT code allowing for all S SPs to be red- 
dened by the same amount of dust described by the Cardelli et al. 
(1989]) law. In our analysis, the spectra and the basis were nor- 
malized to the mean flux in the region 45 80-4620 A. Thus the 
population vector is the fractional contribution Xj of the differ- 
ent SSP models at -4600 A. In addition to the population vector 
the code also returns the fractional contributions fij of each SSP 
to the total stellar mass of the galaxy, which is the more relevant 
physical quantity. The code also returns the velocity shift and the 
Gaussian broadening that need to be applied to the model in or- 
der to fit the observed spectrum. The shifts provide the velocity 
maps for the stars and the broadening is related to the veloc- 
ity dispersion of the stars. From the population vectors we can 
compute the ma ss- and light- w eighted mean age and metallicity 
following |Cid Fernandes et al.| (|2005|) : 



<logf*>L/M = Wj log tj 



7=1 



(Z*)l/m - ^ 



W J Z P 



(2) 



(3) 



where tj and Zj are the age and the metallicity of the j-th SSP 
model, and Wj = Xj or Wj = [ij for light- and mass-weighted 
quantities, respectively. 



3.2.5. Compressed population vectors 



The simulations performed by |Cid Fernand es et al.| ( |2005 



demonstrated that the individual components of the population 
vectors computed by STARLIGHT are very uncertain. Instead 
of analyzing the individual components, |Cid Fernandes et aL 
( 2005 ) showed that a coarsely binned version of the population 
vectors provides a more robust description of the current stel- 
lar c ontent of the galaxies. Thus , follow ing [Cid Fernan des et al. 
(2004]) and |Cid Fernandes et al.| ( 2005] ), we formed compressed 
the population vectors in three age bins: young (age < 300 Myr), 
intermediate (300 Myr < age < 2.4 Gyr), and old (age > 2.4 Gyr) 
stellar populations. 



4. Results 

4.1. Ionized gas 



Figures [2][6] show the main results obtained from the analysis of 
the emission line fluxes. The 2D maps of the galaxy properties 
that are discussed in this section are shown by galaxy in Figs[3][5] 

4.1.1. BPT diagnostic diagram 



Baldwin et al.| ( |1981| ) introduced several diagnostic diagrams 
to segregate spectra of emission-line galaxies and AGNs ac- 
cording to their main excitation mechanism. These diagrams 
are based on easily measured optical emission line flux ratios. 
Figure[2] shows the positions of the galaxies in our sample on 



6 



Stanishev, V. et al.: Type la Supernova host galaxies as seen with IFU spectroscopy 
Table 3. Total galaxy SFR, SFR surface density Xsfr and gas extinction Ay at the SN location derived from our observations. 



SN 


Host galaxy 




total SFR [Mo yr 1 ] 




E SFR at SN position 


A Y at SN position 






Ha a 


STARLIGHT* 


Neill et al.* 


[Moyr^kpc" 2 ] 


[mag] 








<50 Myr 


<0.5 Gyr 








1999dq 


NGC 0976 


5.30 (0.06) 


6.8 


16.0 


10.4 (2.5,36.6) 


8.9(1.2)xl0" 2 


0.97 (0.18) 


200 lfe 


UGC5129 


0.76 (0.02) 


1.0 


3.6 


3.2 (0.6,10.4) 


1.5(0.4)xl0" 2 


0.79 (0.35) 


2006te 


CGCG 207-042 


2.03 (0.05) 


1.6 


5.2 


2.9(1.3,6.9) 


8.3(2.2)xl0" 3 


0.74 (0.33) 


2007A 


NGC 105 NED02 


3.42 (0.04) 


5.3 


9.0 


11.6(2.4,27.5) 


2.6(0.4)xl0" 2 


0.52 (0.22) 


1997cw 


NGC 105 NED02 


3.42 (0.04) 


5.3 


9.0 


11.6(2.4,27.5) 


1.6(0.3)xl0" 2 


0.06 (0.22) 


2007R 


UGC 4008 NED01 


5.33 (0.12) 


5.8 


24.1 


3.6 (2.1,205.6) 


3.2(0.8)xl0" 2 


1.17(0.32) 



(a) this work; from Neill et al. (2009 1. The errors are asymmetric and the values in the parentheses are the + lcr uncertainties. 



the log([Nn] A6584/Ha) - log([Om] A5OO1/H0) diagnostic di- 
agram. The filled circles, filled triangles, and crosses show the 
measurements at the position of the SN, the galaxy nucleus, and 
the total galaxy spectrum, respectively. The dotted and dashed 
lines show two widely used criteria to separate emission-line 
galaxies and AGNs in troduced by jKewley et al.l ([2001 ) and 
Kauffmann et al. ([2003), respectively. From Fig. [2] it is evident 
that the hosts of SNe 2001fe and 2007A/1997cw harbor AGNs 
and the host of SN 199dq is on the border of composite galaxies 
and AGNs. However, the line ratios measured in the total spec- 
tra of these three galaxies still fall into the star-forming region of 
the BPT diagram, which suggests that the AGNs are not strong 
enough to significantly affect the total galaxy spectra. At high 
redshift good S/N, spatially resolved spectroscopy is difficult to 
obtain and weak AGNs may remain unrecognized in slit spec- 
troscopy because typically the whole galaxy falls into the slit. 
The presence of AGNs, even though weak, may still bias the 
metallicity estimation from integrated galaxy spectra. 

The locations of the SNe fall into the (small) region of the 
BPT diagram with the highest density of SDSS galaxies. This is 
the r egion where the metal-rich galaxies are typically found (see, 
e.g.,[Od Ferna ndes et al.|2007] ). High metallicity in this region 
of the BPT diagram is also expected from the Q3N2 method fo r 
metallicity estimation ( |Alloin et al.|1979[|Pettini & Pagel|2004| ). 
These are indications that for the emission line galaxies in our 
sample, the SNe likely exploded in metal-rich environments. 



4.1 .2. Ha velocity field 

The gas velocity maps derived from the Ha emission line show 
smooth gradients and no apparent irregu larities. The analysis 
with the method of Kraj novic et al. ( 2006| ) shows that the veloc- 
ity fields of all five galaxies are consistent with pure disk rota- 
tion without signs of significant disturbances. The cr-maps (not 
shown here) derived form the width of Ha emission line also 
show a simple structure with a single peak at the center. These 
results suggest that the galaxies in our sample are likely relaxed 
systems. 



4.1 .3. Ha flux, extinction, and star formation rate 

On small scales the Ha flux distribution follows the spiral arms 
visible in the broad-band images. Overall, the Ha flux increases 
toward the centers of the galaxies. It can be seen that all six 
SNe in these galaxies are projected onto regions with strong Ha 
emission with fluxes above the galaxy-average. In the late-type 
spirals NGC 976, UGC 5129, and NGC 105 NED02 there is a 
clear Ha flux deficit in the bulge. Interestingly, these are also the 
galaxies that have AGNs. However, these two phenomena are 
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Fig. 2. BPT diagram ( [Baldwin et aT||l981 



The contours show 

the density of SD SS emission line gal axies. The dotted and the 
dashed lines of |Kewley et aE] ( [2Q0TT ) and |Kauffmann et al.|p003> , 
respectively, separate star-forming galaxies, AGNs, and compos- 
ite galaxies. The filled circles, filled triangles, and crosses show 
the measurements at the position of the SN, the galaxy nucleus, 
and the total galaxy spectrum, respectively. 



probably unrelated. Jam es et al. ([2009) studied the radial distri- 
bution of Ha emission in a large sample of spiral galaxies and 
found that the late-type spirals (Sc+) show a Ha flux deficit in 
the bulge regardless of the presence of bars. This effect is much 
less pronounced in the Sa-type spirals or even absent in their 
barred counterparts, which tend to have a high concentration of 
the Ha emission toward the center. This is clearly the case for 
UGC 4008 NED01, which is the only Sa emission line galaxy in 
our sample. 

According to Jam es et aL] ([2009), there is a significant differ- 
ence in the Ha radial profile of barred and unbarred Sb galax- 
ies. The unbarred Sb galaxies show a smooth profile similar to 
Sa galaxies. The barred counterparts have a strong peak of Ha 
emission at the centers, followed by a decrease of the Ha flux, 
before it increases again because of the Ha emission ring at the 
outer radius of the bar. The barred Sb galaxy CGCG 207-042 
shows exactly the same characteristics with a clear Ha emission 
ring at the outer radius of the bar. Thus, the radial Ha emission 
profiles in o ur galaxies are consistent with the findings of [Tames 
|etaLl ( [2009l ). 
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Fig. 3. Upper row: From left to right, the color SDSS image of NGC 976, the observed Ha flux and velocity maps. Lower row: 
The ioniz ation parameter log(L Q, the visual extinction Ay estimated from the B aimer decrement and the metallicity map derived 
with the Pettini & Pagel (2004) 03N2 method. In all maps presented in this paper, x marks the galaxy center and + the SN position. 
The four contour levels overplotted on the extinction and metallicity maps are derived from the Ha map. The four levels are 0.8, 
0.6, 0.4 and 0.2 of the maximum Ha flux. The x,y coordinates are in arcsec with respect to the map centers. The orientation of the 
images is north - up, east - left. 



The gas extinction maps presented in Figs [3][5] also show in- 
crease of the extinction toward the galaxy center. Although the 
extinction maps do not show small-scale structures as clearly as 
in the Ha flux maps, there is a general trend that the extinction 
increases with the Ha flux. This is expected because an increased 
amount of dust is typically observed in the regions of active star 
formation. The total extinction along the lines of sights of the SN 
position is low, except for SN 1997cw (marked with the leftmost 
of the three signs). The extinction along the SN lines-of-sight is 
given in Table[3] 

Table|3] lists the total on-going SFR and the SFR surface 
density, £ S fr, at me positions derived from the extinction- 
corrected Ha flux map. The Esfr values at the SN position are 
consistent with the disk-averaged values for normal spiral galax- 
ies (see F ig. 5 in| Kennicutt 1998). Our values fall in the upper 
half of the Kennicutt (1998) distribution, which can be attributed 
to the SN being projected on regions with higher- than- average 
Ha flux. 

For comparison the total SFR and its confidence intervals 
derived by Neill et al. ( 2009]) ar e also given in Table |3| In general, 



our values are consistent with Neill et al. 



2009 ), although in 



all cases but one we derive lower values. However, it should be 
noted that the values of |Neill et al.| ( |2009] ) were derived with a 
completely different technique - fitting model galaxy SEDs to 
broad-band photometry - and represent the average SFR during 
the last 0.5 Gyr, while our estimates from the Ha flux represent 
the very resent, < 20 Myr, SFR. 



4.1 .4. Ionization parameter and electron density 

The ionization parameter log(£/) - the ratio of the ionizing pho- 
ton density to the gas density - is a measure of the degree of 
ionization of the nebula and can be determined from the ratio 
of two lines of the same element corresponding to two different 
ionization states. The ionization maps of the galaxies were com- 
puted from the rati o of the [Q n] /137 27 and [Om]/15007 lines 
using the relation of Dia z et al. (2000 ). 

The three AGN galaxies clearly show an increased degree 
of ionization toward the center, while the remaining two galax- 
ies do not. In three of the galaxies, NGC 976, CGCG 207-042, 
and NGC 105 NED02, there is also a hint of increasing of the 
ionization parameter toward the outer spirals. The mean ioniza- 
tion parameters for all five galaxies fall into a rather narrow in- 
terval of log(t/) = -3.6 -T- -3.4. The [On] ,13727/ [Om] ,15007 
line ratio is known to provide lower values for the ionization pa- 
rameters compared to other available methods (e.g.,|Diaz et al. 



2000). In comparison with the ionization parameter maps that 
were computed as part of the |Kobulnicky & Kewley| ( [2004 



method for oxygen abundance determinations, the [O n] A3121/ 
[O m] /15007-based maps show very similar features, but are 
shifted toward lower values by ~ 0.3 - 0.4 dex. Even taking 
this offset into account, the average values for our galaxies fall 
int o the l ower e nd of the distribution of the H n galaxies studied 
by |Diaz| ( [l998| ). This part of the distribution is mostly populated 
with Hn galaxies without measurable [Om] /14363 line, which 
tend to be metal-rich. 
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Fig. 4. Same as Fig.^but for UGC 5129 and CGCG 207-042. 



Given the spectral resolution and wavelength range of our & Ferland 2006 ). The ratio of these two lines is sensitive to n Q 
spectroscopy, the electron density, n Q , can be estimated only in the range ~ 10 2 - 10 4 cm -3 . Unfortunately, for the two high- 
from the flux ratio of the [S n] ^6716/6731 lines ( Osterbrock est redshift galaxies in our sample, CGCG 207-042 and UGC 
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Fig. 5. Same as Fig.^but for NGC 105 NED02 and UGC 4008 NED01. 



4008 NED01, the [S n] lines are outside the covered wavelength -1.4 indicate low electron density n e < 10 2 cm -3 and are similar 
range. For the remaining three galaxies the ratio is constant ~1.4 to the measurements of n Q in other galaxies, e.g. the sample of 
across the galaxy and no apparent structure is visible. Ratios of 
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SDSS galaxies studied by |Kniazev et al. (2004). The low elec- 
tron density suggests that there are no shocks in the galaxy. 



4.1.5. ISM oxygen abundance 



Figures [3p] also show the distributions of the metallicity (indi- 
catedas 12+log(0/H) P po4) estimated with the 03N2 method of 
Pettini & Pagel ( 2004 ). For the three galaxies that likely harbor 
AGNs, namely, NGC 976, UGC 5129, and NGC 105 NED02, 
the central regions affected by the AGN are masked. The de- 
crease of the S/N in the outer parts of the galaxies affects 
the measurements of the line ratios and some spaxels are also 
marked as affected by AGN according to the Ke wley etaL] ( |200 1 \ 
criterion. Because no AGN activity is expected in the outer parts 
of the galaxies the line ratios measured in those spaxels are likely 
dominated by the noise and have also been masked in the plots. 
In these three galaxies there are indications for ring-like struc- 
tures with enhanced metallicity and Ha flux. However, it is diffi- 
cult to assess whether these are real structures or artifacts caused 
by the central AGN altering the line ratios. Thus, these ring-like 
structures should be regarded with caution. 

Figure [6] shows the dependence of the metallicity on the 
de-projected galactocentric distance, which was computed from 
the position angle and inclination derived from the analysis 
of the Ha velocity fields. Excluding the host of SN 2007R, 
which has a quite uniform metallicity distribution, the other 
four galaxies show decreasing of the metallicity with the ra- 
dius. The plots also suggest that the high metallicities measured 
in some of the outermost spaxels in UGC 5129, CGCG 207- 
042, and NGC 105 NED02 are most likely due to noise in the 
line flux measurements. The solid lines are guides to the eye, 
showing metallicity gradients of -0.022, -0.035,-0.030 and 
-0.058 dexkpc" 1 for NGC 976, CGCG 207-042, UGC 5129, 
and NGC 105 NED02, respecively. These metallicity gradients 
fall well into the range observed in other nearby galaxies, e.g. 
M51 JBresolin et al.|2004|), NGC 300 ( |Bresolin et al.|2009| ), and 
NGC 628 ( [Rosales-Ortega et al.(2QTT] ). In contrast, in UGC 4008 
NED01 the metatllicity is nearly uniformly distributed with a 
hint of a very small positive gradient of +0.007 dexkpc -1 . The 
blue squares in Fig. [6] show the metalliciy estimates from the az- 
imutally averaged spectra described in Sect. [3] These estimates 
trace the measurements on the individual spaxel spectra very 
well. 

Figure [6] and the inspection of the 2D maps reveals that the 
SN explosion sites are projected onto regions that have the high- 
est, or close to the highest, metallicity within the corresponding 
galaxy. Table[4] shows the metallicity measurements in the total 
galaxy spectra, at the nucleus, and at the SN position. The metal- 
licities at the SN positions in all five galaxies are very similar to 
each other, 12+log(O/H)p P0 4 ~ 8.8 and 12+log(O/H) KK 04 ~ 9.1, 
and are on average by 0. 1 dex higher than the metallicities mea- 
sured from the total galaxy spectra. For the three galaxies that 
host AGNs we also computed the metallicity in the total spectra, 
excluding the central spaxels, which are affected by the AGNs. 
The metallicity was in all cases identical to the one measured 
from the total spectra that included the central spaxels (the lat- 
ter are given in Table]?]). This result shows that in these cases 
the AGNs are too weak to significantly affect the total galaxy 
spectrum and the metallicity estimation. 
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Fig. 6. Radial dependence of 12+log(0/H)p P o4. The red dots 
show the measurement at the SN position. The horizontal dashed 
lines indicate the metallicity measured from the total galaxy 
spectra (Table [4]). As a guide to the eye, we plot metallicity gra- 
dients of -0.022, -0.035,-0.030, -0.058 and +0.007 dexkpc" 1 
(from top to bottom) with the solid lines. We estimate that the ac- 
curacy of these gradients is ~ 0.005. The blue squares show the 
metallicities measured from the azimuthally averaged spectra. 
The points shown with the cross and plus signs are the metallic- 
ities of the H ii regions discussed in Sect. |4.1.6 
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Table 4. ISM metallicity estimates from the total galaxy spectra, at the nucleus, and at the position of the SN using the PP04 and 
KK04 methods. 



SN 


Host galaxy 






12+log(0/H) 






total galaxy spectrum 


nucleus 


at the SN position 


PP04 


KK04 


PP04 


KK04 


PP04 


KK04 


1999dq 


NGC 976 


8.75 (0.17) 


9.08 (0.05) 


8.69* (0.08) 


8.84 a (0.23) 


8.81 (0.05) 


9.15 (0.03) 


2001fe^ 


UGC5129 


8.72 (0.16) 


8.98 (0.11) 






8.82 (0.10) 


9.10(0.09) 


2006te 


CGCG 207-042 


8.57 (0.14) 


8.96 (0.10) 


8.72 (0.06) 


9.10(0.05) 


8.75 (0.09) 


9.03 (0.13) 


2007A^ 


NGC 105 NED02 


8.65 (0.14) 


9.02 (0.08) 






8.83 (0.06) 


9.09 (0.07) 


1997cv/ 


NGC 105 NED02 


8.65 (0.14) 


9.02 (0.08) 






8.76 (0.06) 


9.10(0.04) 


2007R 


UGC 4008 NED01 


8.75 (0.16) 


9.08 (0.06) 


8.73 (0.12) 


8.66 (0.55) 


8.73 (0.09) 


9.12(0.06) 



(a) 



may have AGN contamination. (b) the metallicity not measured at the center because these galaxies harbor AGNs. 



4.1 .6. Notes on the individual galaxies 



NGC 976 : The metallicity distribution is nearly symmetric 
around the galaxy center except for a somewhat extended region 
located at coordinates (+8, +22). The metallicity of this region is 
higher compared to the other parts of the galaxy at the same ra- 
dial distance. Examining the [Om]5007/H/3 and [Nll]6584/Ha 
maps (not shown in the paper), we noted that this is caused by an 
asymmetry in [O lll]5007/H/3. Both ratios are nearly symmetri- 
cally distributed around the galaxy center except for the region 
at (+8,+22), which has a lower [O iii]5OO7/H0 ratio resulting in 
a higher metallicity estimate. There is also a slight decrease of 
the degree of ionization at the same location. The galaxy was in- 



cluded by Maio lino et al. ( 1997| ) in the control sample for their 



study of the cause of the elevated star formation in Seyfert 2 
compared with Seyfert 1 galaxies. The authors found no obvi- 
ous trigger of star formation in NGC 976. 
NGC 495: This red barred Sa galaxy shows no emission lines. 
Mill er et al. (2002) found it to be a member of a poor galaxy 
cluster, which was the richest cluster among those studied in 
their work, however. It is therefore possible that the gas com- 
ponent of NGC 495 was separated from the galaxy by the tidal 
interaction with the other cluster members. 
UGC 5129: Thi s galaxy was includ ed in the study of isolated 
disk galaxies by |Varela et al.| ([2004]). It was included in the fi- 



nal list of 203 galaxies (out of an initial 1706) that were likely 
not affected by other galaxies during the last few Gyr of their 
evolution. 

CGCG 207-042: The spirals arms of the galaxy are barely vis- 
ible in the SDSS image. However, there are three Hn regions 
along one of them that are clearly visible in the Ha map. They 
are roughly located at (x, y) coordinates (+8,-27), (+27,-10) and 
(+20, +27). The metallicity decreases considerably along the spi- 
ral arm, which is also accompanied by a strong increase of the 
ionization parameter. The two outermost H n regions are shown 
with blue and magenta points in Fig. [6] 

NGC 105 NED02: In the Ha velocity map there is a spot located 
at (+28, +4) that clearly does not follow the velocity of the un- 
derlying part of the galaxy but moves away -120 km s" 1 faster. 
At the same position there is a very faint spot in the broad-band 
images. This sport also clearly shows increased Ha emission and 
a marginal increase of the ionization parameter. The metallicity 
of the spot is lower than the rest of the galaxy by at least 0.2 dex. 
The points corresponding to this spot are shown with blue points 
in Fig. [6] Given the properties of this feature, it is possible that 
this is a dwarf satellite galaxy of NGC 105 NED02. Another in- 
teresting feature is that the central ring-like pattern of increased 
metallicity is interrupted by a region of slightly lower metallic- 
ity located at coordinates (-3,-2). Examining the [O m]5007/H/3 



and [Nll]6584/Hof maps, we again noted that this is caused by 
an asymmetry in [O iii]5OO7/H0. The [O m]5OO7/H0 ratio at (- 
3,-2) is slightly higher and causes the lower metallicity estimate. 
UGC 4008 NED01: This is the only galaxy in our sample that 
shows a positive metallicity gradient. 

4.2. Stellar populations 

4.2.1 . Stellar vs. gas dust extinction 

STARLIGHT fits provide an estimate of the extinction by dust suf- 
fered by the stellar light. The assumption that the stellar popula- 
tions of different age are subject to the same extinction is proba- 
bly not entirely correct. It is reasonable to assume that the young 
populations can be still embedded in the dusty nebula where the 
stars formed and can be subject to higher extinction. STARLIGHT 
has the capability to take this into account and can determine dif- 
ferent extinctions for the different SSP models. However, this ap- 
proach adds additional uncertainty to the already complex prob- 
lem of recovering the properties of the stellar populations. We 
have chosen to assume a single extinction for all SSPs. 

The extinction maps of the star light are shown in Fig. [7] and 
there were no easily identifiable features in them. In compari- 
son with the extinction derived from the emission lines (Figs [3} 
[5]), the extinction derived by STARLIGHT fits is lower. The re- 
lation between the star and gas extinction shows considerable 
scatter and the two quantities appear to be uncorrelated, except 
for NGC 105 NED02. In this galaxy there is a clear linear rela- 
tion between the star and gas extinction, with the gas extinction 
being about tw ice the star's extinct i on. A similar relation was 
also derived by |Cid Fernandes et al.| ( |2005] ) in their analysis of a 
sample of SDSS galaxies. 

4.2.2. Mean stellar age and metallicity 

Table[5] lists the mass- and light- weighted mean stellar popula- 
tion age and metallicity determined from fitting the total galaxy 
spectrum formed as the sum of all spaxels with (these fits are 
shown in Fig. [5]) and without the AGN-affected central spax- 
els. The results show that all galaxies in our sample have a 
higher mean stellar metallicity than solar. This is in accord with 
the findings from the emission lines analysis. The mean mass- 
weighted stellar age of the five emission line galaxies is ~ 5 Gyr. 
NGC 495, which shows no emission lines, has an older stellar 
population of about 12 Gyr. These values can be compared with 
studies based on total galaxy spectra, e.g. obtained w ith drift- 
scanning with a long-slit of local galaxies (e.g., [Gallagher et al.| 
2005 ) or spectroscopy of high-redshift galaxies when practically 
the whole galaxy light falls into the slit. Note that the residuals 
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Table 5. Stellar population metallicity and age estimates from the STARLIGHT fits. 



SN Host galaxy total galaxy/minus AGN a at the SN position^ 







<loga*)) M 


<z*> M 


<loga*)> L 


<z*> L 


<logfe)> M 




<log(*.)> L 


(Z*) L 


1999dq 


NGC 976 


9.89/9.96 


0.039/0.036 


8.82/8.80 


0.029/0.030 


9.75 


0.042 


8.82 


0.030 


1999ej 


NGC 495 


10.17 


0.045 


9.96 


0.038 


10.14 


0.042 


10.03 


0.037 


2001fe 


UGC5129 


9.66/9.96 


0.038/0.023 


8.95/8.92 


0.029/0.026 


9.45 


0.037 


8.71 


0.031 


2006te 


CGCG 207-042 


9.46 


0.029 


8.56 


0.022 


9.55 


0.027 


8.83 


0.022 


2007A 


NGC 105 NED02 


9.86/9.80 


0.030/0.033 


8.53/8.51 


0.026/0.031 


9.78 


0.040 


8.48 


0.029 


1997cw 


NGC 105 NED02 


9.86/9.80 


0.030/0.033 


8.53/8.51 


0.026/0.031 


9.77 


0.039 


8.43 


0.028 


2007R 


UGC 4008 NED01 


9.68 


0.042 


9.01 


0.031 


9.96 


0.044 


9.11 


0.030 



(a) 'total galaxy' - values derived by fitting the spectra formed by summing (un- weighted) all spectra in the data cubes, 'minus AGN' - with the 
AGN affected spaxels excluded. w interpolated from the maps in Fig.[7] 



show a large-scale pattern with a full amplitude of up to -4% 
(Fig. [8]). This signals either a problem with the relative flux cali- 
bration of the observed spectra or a problem in the SSP models. 
At present it is difficult to quantify what effect this would have 
on the results that are based on the spectral fitting. 

Figure [7] shows the mass- and light-weighted mean stel- 
lar population age and metallicity maps of the six galaxies. In 
Table[5] are given the measurements for the total galaxy and at 
the SN position. The maps show considerable scatter and it is 
difficult to identify clear structures in them. Many spaxels that 
indicate high metallicity appear in the outer parts. This is most 
likely not real but rather a result of the insufficient S/N of the 
spectra even after applying the Voronoi binning. Nevertheless, 
there may be a slight increase of the metallicity toward the cen- 
ter, especially in NGC 495 and UGC 4008 NED01. The same is 
also true for the stellar age maps, and again there is an indication 
of an older stellar population toward the nucleus, which can be 
expected. 

To investigate the matter in more detail, we plot in Fig. [9] 
the the mass- and light-weighted mean stellar population age 
and metallicity measurements as a function of the de-projected 
galactocentric distance. Unfortunately, the plot confirms that the 
measurements from the individual spaxels spectra show too large 
scatter. Unlike the ionized gas metallicity measurements, which 
show a small scatter of <0.05 dex at a given radius (Fig. [6]), 
the stellar metallicities estimated from the STARLIGHT fits show 
scatter as large as 0.3 dex, for example at a radial distance of 6 
kpc in NGC 976, UGC 5129, and NGC 105 NED02. The age 
estimates also show considerable scatter. 

The analysis of the emission lines shows that most of the 
ISM properties have a well-defined axial symmetry. One can ex- 
pect this to be also the case for the stellar populations and hence 
asymmetries are unlikely to be responsible for the observed scat- 
ter in the outer parts of the galaxies. The scatter clearly increases 
with the radial distance (Fig. [9}, suggesting that the lower S/N 
of the spectra in the outer parts of the galaxies is causing it. 
After the Voronoi binning the analyzed spectra a have minimum 
S/N~ 15-20 at 4600A. The large scatter that we observe in the de- 
rived quantities demonstrates the limitations of the full- spectrum 
fitting technique in the low- S/N regime and suggests that an S/N 
significantly higher than 20 is needed to achieve reliable results. 

Given the large scatter of the measurements from the indi- 
vidual spaxel spectra, interpolating at the location of the SNe 
from the 2D maps is not recommended. An alternative approach 
is to use the measurements obtained from the azimuthally aver- 
aged spectra and interpolate them at the radial location of the 
SNe. This approach is better when there is evidence that the 
galaxy properties are symmetric around the nucleus. In Fig. [9] 
the blue symbols show the values estimated from the fits of the 



azimuthally averaged spectra and the vertical dashed lines show 
the radial distance of the SNe. The corresponding fits are shown 
in Figs. [T0|[T2l The mean age and metallicity show a smooth 
radial dependence. In some cases the metallicity derived from 
the azimutally averaged spectra suggests negative gradients of 
up to -0.03 dex kpc -1 . However, given the large uncertainty 
with which the stellar metallicity is estimated (> 0.2 dex), the 
significance of these gradients is difficult to assess. The mean 
ages qualitatively show the same behavior with decreasing age 
outward. We note that the different types of weighting, mass 
or light, lead to different radial dependencies, with the light- 
weighted quantities showing stronger variation. The metallic- 
ity and the age at the locations of the SNe linearly interpolated 
from these radial dependencies are given in Table[5] It is also 
worth mentioning that the light-weight quantities appear to have 
a slightly lower scatter, most pronounced in the inner regions 
where the spectra have a higher S/N. In most cases the light- 
weighted metallicies are lower than the mass- weight ones. The 
light-weighting gives much more weight to the younger stellar 
population and this result may imply that the younger popula- 
tions have lower metallicity. 

4.2.3. Binned population vectors 

The 2D maps representing the fractional contribution of the 
young (age < 300 Myr), intermediate (300 Myr < age < 2.4 Gyr), 
and old (age > 2.4 Gyr) stellar populations are shown in Fig.[T3| 
In Fig. [14] the measurement from the individual spaxel spec- 
tra and the azimutally averaged spectra are plotted vs. the de- 
projected galactocentric distance. The measurements from the 
individual spectra again show considerable scatter in the outer 
parts of the galaxies. For this reason we again estimated the val- 
ues at the radial distance SN from the azimutally averaged spec- 
tra and not from interpolation of the 2D maps. The estimated 
stellar population fractions at the SN radial distances are given in 
Table[6] along with the values derived from the total galaxy spec- 
tra. The S/N of the spectra used to derive these values are also 
shown. The analysis shows that the five emission line galaxies 
contain stellar populations of different ages, including a consid- 
erable fraction of young stars. In general, there is a clear trend 
of increasing the fraction of young stars with the radial distance. 
Depending on the galaxy, at a distance of 4-8 kpc the trend is re- 
versed and the fraction of young stars starts to decrease. In four 
of the galaxies the fraction of old stellar populations monotoni- 
cally increases toward the galaxy nucleus, which is expected for 
most star- forming spiral galaxies. The exception is CGCG 207- 
042, the host of SN 2006te, which shows a decrease of the frac- 
tion of old stellar populations toward the center. NGC 495 is 
dominated by old stellar populations with possibly a small frac- 
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Table 6. Compressed population vectors showing the contribution of the young (age < 300 Myr), intermediate (300 Myr < age < 
2.4 Gyr), and old (age > 2.4 Gyr) stellar populations to the formation of the observed total galaxy spectrum. The total galaxy stellar 
masses derived from our STARLIGHT fits and by |Neill et al.| ( |2009] ) are also given. 



SN 


Host galaxy 




total spectrum 




at the SN position^ 7 




log(M* [M Q ]) 






Young 


Inter. 


Old 


S/N fl Young 


Inter. 


Old 


S/N c 


this work 


Neill et al. (2009 ) 


1999dq 


NGC 976 


0.19 


0.49 


0.32 


81 


0.19 


0.41 


0.40 


84 


10.98 


10.78 


1999ej 


NGC 495 


0.00 


0.11 


0.89 


54 


0.00 


0.00 


1.00 


33 


10.85 




2001fe 


UGC5129 


0.16 


0.40 


0.44 


97 


0.29 


0.33 


0.38 


81 


10.22 


10.22 


2006te 


CGCG 207-042 


0.25 


0.35 


0.40 


35 


0.21 


0.18 


0.62 


42 


10.25 


10.31 


2007A 


NGC 105 NED02 


0.31 


0.35 


0.34 


78 


0.36 


0.38 


0.26 


82 


10.61 


10.87 


1997cw 


NGC 105 NED02 


0.31 


0.35 


0.34 


78 


0.38 


0.38 


0.24 


82 


10.61 


10.87 


2007R 


UGC4008NED01 0.12 


0.36 


0.53 


54 


0.14 


0.29 


0.57 


85 


11.10 


10.98 


(«) s/n of the total galaxy spectra. {b) values at the SN radial distance 


interpolated from the radial depend 


encies derived from the azimul 


averaged spectra (Fig.[l4j). (c) 


S/N of the azimuthally avera| 


*ed spectrum that is closest to the radial distance on the SN. 




tion of younger stars in the central few kpc. The contribution of 


4.2.5 


Current stellar mass 







the younger population is small, however, and its presence can- 
not be confidently confirmed. The behavior of the intermediate 
age stellar populations is the opposite to that of the old ones. 



Cid Fernandes et al. ( 2005 ) showed that the compressed pop- 



ulation vectors can be recovered with an accuracy better than 
~10% for S/N>10. However, considering the uncertainties in- 
volved in the computation of the SSP models as well as other 
uncertainties such as the correlations between the fitted parame- 
ters, the relative flux calibration and the dust extinction laws in 
the galaxies, the accuracy is probably no better than -10%. This 



is also supported by the level of the scatter in Fig. 14 In this 
context, the population vectors at the locations of SNe 1999dq, 
2007 A, 1997cw, and 2007R are un-distinguishable from those 
of the whole galaxies (Table[6]). At the location of SN 2006te 
there is a larger contribution from old populations at the expense 
of the intermediate age, while the fraction of young stars is the 
same as for the whole galaxy. For SN 200 lfe there is marginal 
evidence for an increased contribution of a young population at 
the position of the SN. The host of SN 1999ej formed the bulk of 
its stars about 13 Gyr ago followed by a less intense star-forming 
period about 2 Gyr ago. At the distance of SN 1999ej we only 
find evidence for the older population. 



4.2.4. Stellar kinematics 

The velocity dispersion maps derived from the STARLIGHT fits 
show a simple morphology with a single peak centered at the 
galaxy nucleus. In Fig. 15 the velocity fields of the stars in the 
five emission line galaxies are compared to the velocity fields 
derived from the Ha emission line. The two maps are very simi- 
lar and small systematic differences are only revealed after sub- 
tracting the two maps (the last column in Fig. [T5] ). Evidently, 
the gas rotates faster in the central regions than the stars, with 
the difference being largest in UGC 4008 NED01. These dif- 
ferences between the rotation of stars and ionizied gas in the 
central regions of galaxies are well-known and have been ex- 
tensively studied (see, e.g., Pizzella et al.|2004 and references 



therein). We note that none of the galaxies shows a sign of 
coun ter-rotating gaseous disk (e.g., Ber tola et al.||1996| |Rubin 
|et al-| 1992| ). 

The stellar velocity fields were also analyzed with the meth- 
ods of Krajnovic et al. ( 2006| ). As with the Ha velocity map, 
within the uncertainty we also found no evidence for deviations 
from pure disk rotation. 



An estimate of the present-day stellar mass of the galaxies was 
obtained from the STARLIGHT fits of the total galaxy spectra. 
The fits are shown in Fig. [8] and the masses are given in Table 6] 
All the galaxies have masses exceeding 2 x 10 10 Mq and can be 
classified as quite massive. The fairly high metallicity that we 
derived for both the ionized gas and the stellar component are 
therefore in line with the expectation from the mass-metallicity 
relation, e.g. Tremonti et al. ( 2004| ). We note that the va lues that 
we obtain are very close to those of Neill et al.| ([2009]), which 
were obtained by a different methodo^ 



ogy (see Sec. |4. 1.3) . 



5. Discussion 

5.1. Galaxy mass and metallicity 

We used IFU spectroscopy to derive the spatially resolved prop- 
erties of six face-on spiral galaxies that hosted seven nearby 
SNe la. The masses of the galaxies derived from the analysis 
of the total spectra wit h the STARLIGHT code are all higher than 



2 x 10 1Q M G . Recently, |Kelly et al.| ( |2QlQ| ), |Sullivan et al.| ( |2Q10| , 
and |Lampeitl et al.| ( |2010| ) have claimed that the residuals from 
the best-fit Hubble line corr elate w ith the SN host stellar mass. 
Furthermore, Sulliv anet al.| ([2010) proposed to incorporate into 
the cosmological SN la analyses two different absolute peak 
magnitudes for SNe in hosts with masses lower or higher than 
10 10 Mq; after the "lightcurve width - luminosity" and color 
corrections the SNe in the more massive hosts are found to be 
~ 0.06 - 0.09 mag brighter than their counterparts in lower mass 
hosts. The galaxies in our sample fall into the high-mass/low- 
specific SFR bins defined by Sulli van et al. ( 2010 ). Accordingly, 
one can expect the SN in these galaxies to have on average neg- 
ative Hubble residuals. From Table |2] one can see that four of 
the SNe have significant positive residuals (> 2<x). The other 
three have negative residuals, but only one of them is bigger 
than the uncertainly. The mean weighted residual is positive, 
+0.07 ± 0.22; however, one should keep in mind that we used 
only very few SNe in our analysis. 

The cause of the apparent dependence of the SN la lumi- 
nosity on the host galaxy stellar mass is still unclear. Theoretical 
investigations have shown that various parameters of the explod- 
ing WD, such as its metallicity, C/O ratio, central density, and 
progenitor age can affect the amount of 56 Ni synthesized in the 
expl osion to a different degree and hence the SN luminosity (see, 
e.g.,|Bravo et al.pOTOl powell et al.||2009| |Ropke et al.||2506l 
Timmes et al. 2003, and references therein). Among these pa- 
rameters, the metallicity is known to correlate with the galaxy 
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Fig. 16. Mass-weighted stellar metallicity vs. gas-phase oxygen 
abundance estimated from the azimutally averaged spectra. 



mass (see, e.g., Tremonti et al. 2004) and is likely to have the 
strongest impact. Our analysis of the emission line fluxes and the 
stellar populations revealed that the galaxies in our sample have 
on average solar and higher metallicity (Tables [4] and [5]). This is 
not surprising because the galaxies are quite massive a nd by the 
virtue of the mass-metallicity relation (see, e.g., |Tremonti et al. 



2004) may be expected to have high metallicities. For five of the 
SNe, the ISM metallicity measured at the location of the SN is 
higher than the galaxy average by about ~ 0.1 dex (Table [4]). 
This can be explained by the presence of radial metallicity gra- 
dients and our target selection criteria. Figures [6] and [9] show that 
the galaxies in our sample have radial metallicity gradients. At 
the same time, the selection criterion that the SNe are located 
on a high surface brightness location in the galaxies led to a 
SN sample that is biased toward SNe close to the galaxy nu- 
clei. Together with the presence of the metallicity gradients, this 
resulted in most of the SNe being at locations wit h higher-than- 
average metallicity within the galaxies (see also Kewley et al. 
|2Q05> . 

While the gas-phase metallicity is easier to estimate, a more 
relevant quantity is the stellar metallicity. Figure 16 shows the 
mass- weighted stellar metallicity vs. the gas-phase oxygen abun- 
dance estimated from the azimutally averaged spectra. The two 
quantities appear to be correlated. The slope of the li near fit is 
~ 1.8 with a dispersion of ~ 0.1 dex. Note that Cid Fernandes 
et al.| p005| ) also found that the gas-phase and the stellar metal- 



licities are correlated from an analysis of a large sample of SDSS 
galaxies. However, the two relations are difficult to compare be- 
cause the Cid Ferna ndes et al.| ( [2005] ) analysis also included low- 
metallicity galaxies and galaxies in a somewhat higher redshift 
interval. 



5.2. The impact of the metallicity gradient 

The presence of abundance gradients in both s piral and elliptical 
galaxies is now a well-established fact (e.g., |Henry & Worthey 
1999; Zaritsk yet al.| [T994). If not taken into account, the gra- 
dients will affect any attempt to study the properties of SNe la 
and/or their progenitors as a function of their host galaxy metal- 
licity. The values of the gradients seen in the galaxies in our sam- 
ple suggest that SN la progenitors that form at radial distances 
greater than ~ 15 kpc may have metallicities that are lower by 
a factor at least 2-3 than progenitors in the central parts. Studies 
of the radial distribution of SNe la within their hosts galaxies 



have shown that more SNe explode in the central regions (e.g., 
I Bartunov et al.|[2007] |Ivanov et aL|2000| |van den Bergh][T997l 
Wa ng et al.||1997| )rHowever, SNe la are also found at large 
galactocentric distances in both spiral and elliptical galaxies. In 
Fig. 17 we show the distribution of the projected galactocentric 
distances (PGD) for a sample of 305 SNe with modern CCD 
observations (observed after 1990) and with known host galaxy 
type, redshift, and offset from the center. About 7% of the SNe in 
spiral galaxies and 20% in the ellipticals are found at PGD> 20 
kpc. Since the real galactocentric distances are always greater 
than, or equal to, the PGD, the above-mentioned fractions are 
lower limits. Therefore, a significant fraction of SNe may have 
progenitors with a metallicity that is much lower than that of the 
host average. 

An important question is whether the present-day galaxy 
metallicity is a good proxy of the met allicity of SN la progen - 
itors. This was recently studied by |Bravo & Badenes| p011| ), 
who used simplified one-zone galaxy evolution models coupled 
with the SN del ay-time distribution (DTD) functions of Pritchet 
|et al. | ( 12008] ) and |Maoz et~aL] ( |2010| ). The authors concluded that 
the galaxy ISM metallicity is a good proxy for the SN progen- 
itor metallicity and derived simple linear relations to estimate 
the progenitor metallicity from the present-day host metallicity. 
However, |Bravo & Badenes (2011) did not include metallicity 
gradient, and more importantly, its possible evolution with time. 
There is growing evidence that the disks in late-type galaxies 
formed and evolved slowly under the constant inflow of metal- 
pool gas from the galactic halo. The galaxy chemical evolution 
models and hydrodynamical simulations have shown that the 
metallicity gradient evol ves considerably during the last 10 Gyr 
of the galaxy evolution ([Chiappini et al.||200l[ |Fu et al.|[2009 



Molla & D iaz 2005 ; |Molla et al.||1997t|Naab & Ostrike r|2006 
Pilkingto net al.|2012[ |Sanchez-Blazque z et al.|2009j ). Although 



the exact results d epend of the particular code and model used 
(Pilkington et al. |2012| ), all studies but one ( [Chiappini et al. 
200 1 J) show that the metallicity gradient was steeper in the past 
and gradually flattens out to reach present-day values similar 
to those observed in local spiral galaxies. Recently, there has 
also been obs ervational suppor t for this conclusion. |Yuan et al. 
( 2011 ) and |Jones et al. ( |2010| ) reported metallicy gradients of 
-0.16 dex kpc" 1 and -0.27 dex kpc -1 for galaxies at redshifts 
z=1.5 and z=2.0, respectively. We note that the galaxy chem- 
ical evolution studies show that the mean disk metallicity has 
increased slowly by ~ 0.3 - 0.5 dex during the last several Gyr. 
The gradients seen in the galaxies in our sample and in other 
galaxies at low and high redshift imply that the metallicity dif- 
ferences within the same galaxy may exceed the cosmological 
increase of the mean metallicity. In addition, some studies have 
pointed out that the metallicity gradient in the outermost parts of 
the galaxies may be steeper than in the inner disk (see, e.g., Fu 
et aL|2009> . 

The above studies highlight the complexity of estimating 
the metallicity of the SN la progenitors from their host galaxy 
present-day metallicity. The difference between the present host 
metallicity and the SN progenitor metallicity is a complex func- 
tion of several factors, some of which are poorly understood and 
not very well constrained with observations: the radial distance 
at which the progenitor formed, the age of the progenitor, and 
the evolution of the metallicity gradient. For example, a progen- 
itor that formed at large radial distance will have increasingly 
larger difference from the preset-day metallicity at the same ra- 
dius as the progenitor ages. Another uncertainty can be added 
if the galaxies have experienced major mergers and radial star 
migrations, which tend to flatted the metallicity gradient (see, 
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Fig. 17. Distribution of the observed galactocentric distance for 
a sample of 305 nearby SNe la in late- and early-type hosts. 
Because the distances have not been de-projected these are the 
minimum galactocentric distances. 



e.g., |Kewley et aLl|2010[ |Sanchez-Blazquez et aL]|2Q09| ). This 
complexity may be the reason why the attempts to correlate 
the Hubble residual with the host global metallicity have not 
led to conclusive results (|D T Andrea et al.|2011[ [Gallagher et al. 
200 51 [20081 IHowell et al.|2009) . Note however that |Howell et al. 
(2009} did not directly measure the metallicity but rather esti- 
mated it from a mass-metallicity relation. 

All SNe in our sample but one are within 5 kpc from the 
galaxy centers. Generally, the chemical evolution models show 
that the metallicity close to the galaxy nuclei changes least. 
Therefore, the metallicity of the SN progenitors that formed near 
the center should be closer to the present-day galaxy metallic- 
ity compared to the progenitors that formed in the outer parts. 
Together with the fact that we measured fairly high present-day 
metallicity at the locations of all SNe, this suggests that their pro- 
genitors did not form in metal-poor environments, unless they 
came from very old stellar population with a long delay time. 

5.3. Star formation history 

Much of the recent progress on the question of SN la progenitors 
has been achieved through studies of the SN rates. It is now well- 
established that the SN la rate depends on both th e total stellar 
mass and the recent SFR in the host galaxy (e.g., Brandt et al. 



2010; Mannu ccTet al.||2005[ [Maoz et al.||20T0l |Scannapieco~& 
Bildsten 2005 ; Sulliva net al.|20 06), which led to a two compo- 
nent model for the SN la rate, the co-called A+B model. Along 
with the fact that SNe la are also observed in old, passive galax- 
ies, this points to the existence of at least two evolution channels 
for SNe la associated with young and old stellar populations. 

Except for NGC 495, all other galaxies in our sample con- 
tain a considerable fraction of young stars and strong Ha emis- 
sion, indicating ongoing star formation activity. The STARLIGHT 
fits of the total galaxy spectra are shown in Fig. [8] Except for 
NGC 495, all other galaxies show a similar pattern, namely, the 
population vectors Xj show the largest contribution from SSPs 
with ages 0.5-5 Gyr. Young populations, ~50 Myr, are also con- 
fidently detected in all cases. From Figs. [13] and [9] it can be seen 
that the fraction of young stars increases with increasing the ra- 
dial distance. It is known that the full- spectrum fitting techniques 
tend to estimate sus piciously large components w i th ages ~ 1 
Gyr (see, e.g., |Asari et al.|2007[|Mathis et al.|2006[|Panter et al. 



2007 ) and one may ask whether the large contribution of SSPs of 
similar ages that we see in our analysis is real. |Asari et al. (2007) 
and |Cid Fernandes et a"L] ( |2009| ) report that the problem disap 



peared once they switched from the original |Bruzual & Chariot 
( 2003) fitting basis based on STELIB to a new basis that uses the 
MILES spectral library. Because we also used the newer Bruzual 
& Chariot basis based on MILES, our results are also likely un- 
affected by the above-mentioned problem. 

We estimated the current SFR rate from the Ha emission 
line flux (Table [3]). Because most of the ionizing photons are 
produced by massive, short-lived stars, the Ha flux is a tracer 
of the very recent star formation, < 20 Myr. Another estimate of 
the SFR can be obtained from the STARLIGHT fits following the 
methodology described in Asar i et al. ( |2007| ). We estimated the 
mean SFRs during the last 0.5 Gyr and last 50 Myr. The values 
are given in Table [3] The mean SFRs over the last 50 My r are 
very similar to the estimates obtained from the Ha flux; |Asari| 
|et al.| ( [2007] ) have already demonstrated that there is a tight cor- 
relation between these two estimations using a large sample of 
SDSS galaxies. On the other hand, the mean SFRs during the 
last 0.5 Gyr are by a factor 3-5 higher than the SFR estimates 
from Ha flux, but are similar to those of |Neill et al.| (2009). 
The only exception is UGC 4008 NED01, for which we obtain 
a much higher value; note, however, that the confidence inter- 
val quoted by |Neill et al.| ( |2009] ) has an upper limit higher by an 
order of magnit ude than our estimate. Note also that the model 
SEDs that |Neill et al. 1(2009) fitted to the broad-band photom- 
etry are based on eight galaxy models with pre-defined SFHs, 
which were meant to represent the Hubble galaxy types plus one 
star-burst galaxy model. STARLIGHT does not assume any pre- 
defined SFH and the contributions of all SSPs are free param- 
eters. Therefore, STARLIGHT is much more flexible to describe 
galaxies with arbitrary SFHs. 

If the two components of the A+B model represent the con- 
tribution of two different channels to produce SNe la, we can es- 
timate the probability from which channel the SNe in our galax- 
ies were produced. With the A and B constants estimated by 
Sull ivan et al.| ( |2006] ) and our measurement of the galaxies' total 
mass and SFR, the SNe in the five star-forming galaxies have 
an about equal chance to have come from the young or the old 
channel. 



5.4. Correlation between the SN and host galaxy properties 

Despite the low statistics of our sample, we can still test corre- 
lations between the SALT2 x\ parameter and the Hubble resid- 
uals of the SNe with the various parameters that we derived for 
the total galaxy and at the SN locations (Tables [3][6]). We found 
no statistically significant correlations. The small number of ob- 
jects is certainly not ideal for this analysis, but the metallicity 
and the masses of the galaxies in our sample also span quite a 
small range. To search for correlations it is necessary to expand 
the sample toward lower masses and metallicities. This may not 
be an easy task because SNe la in metal-poor galaxies in the 
local Universe are rare. Besides, the low-mass, low-metallicity 
galaxies tend to be faint and are more difficult to observe with a 
sufficient S/N. 

In addition to the above concerns, when correlating the SN 
properties with the properties of their host galaxies at the loca- 
tion of the SN one should always b ear in mind that t he SN pro- 
genitors may be very old stars (e.g., |Maoz et al.|2()TT ). The pro- 
genitor system may have migrated from its birth place and the 
galaxy properties at its present location may be different from 
those where the progenitor has formed. [Maoz et aL ( 2011| ) dis- 
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cussed that random stellar motions will affect the SN progenitor 
and its surroundings in the same way. As a result, the popula- 
tion that gave birth to the SN progenitor will also be present 
in the new SN location. Following the same argument, the ra- 
dial star migra tions should not have a sign ificant effect either. 
Additionally, Sanchez-Blazq uez et al. ( |2009| ) found that within 
the central 5-10 kpc the radial star migration during the whole 
galaxy evolution is fairly small ~1.5 kpc and increases to only 
- 3.5 kpc in the outer parts. 

Another problem is related to the projection effects. Because 
the SN la progenitors have a broad range of ages between 
-100 Myr and 10 Gyr, SNe can explode anywhere along the 
line of sight. An SN produced by an young progenitor would 
have most likely exploded in the galactic disk, where the newly 
formed stars and the ionized gas typically reside. In this case the 
progenitor star's metallicity should be close to that of the ion- 
ized gas at the (projected) location of the SN. In the case of old 
progenitors, however, the SN may have exploded in the galactic 
halo and the metallicity of its progenitor may be very different 
from the gas. Similarly, the stellar continuum at a given spaxel 
is the sum of all star light along the light of sight. All these ef- 
fects indicate that the correlation of the SN properties with the 
properties of its local environment is not unambiguous. 

Clearly, to study the correlation between the SN properties 
and its local environment, a large, unbiased sample of galax- 
ies observed with a large-field IFU spectrograph is needed. Our 
sample of seven SNe la in six galaxies is not large e nough to 
draw any conclusion. The ongoing CALIF^ survey (Sanchez 
et al.||2012| ) will provide IFU observations of about 600 galax- 
ies at redshift z ~ 0.02. CALIFA uses the same instrument as the 
observations presented in this paper with similar setup and expo- 
sure times, and will provide data of similar quality as ours. The 
CALIFA targets are selected from a larger pool of about 1000 
galaxies based only on the visibility of the targets at the time 
of the observations. Many of these galaxies are known to have 
hosted SNe. In addition, there are several ongoing large-field SN 
searches with that will potentially discover many new SNe in 
CALIFA-targeted galaxies. Thus, the CALIFA survey will pro- 
vide a solid base to further expand the studies SNe la proper- 
ties as a function of their local environment to all SN types. In 
addition, CALIFA will provide the full galaxy spectra that will 
allow one to avoid the aperture effects to which the SDSS spec- 
troscopy is subject. At low redshift the 3 " -diameter fibers of the 
SSDS spectrograph cover only the galaxy nucleus, whose prop- 
erties may be very different from those of the disk and may not 
be representative for the environment of the SNe that exploded 
elsewhere in the galaxy, for example because of the radial gradi- 
ents of these properties. 

6. Conclusions 

In this pilot study we have obtained and analyzed IFU spec- 
troscopy of six nearby spiral galaxies that hosted seven SNe la. 
For the data reduction we developed and tested a robust reduc- 
tion pipeline. A set of tools that implement various methods to 
derive the properties of the ionized gas and the stellar popula- 
tions from the data-cube were also developed. This allowed us 
to generate 2D maps of the galaxies properties. The analysis of 
the maps showed that the quality of the data is sufficient to accu- 
rately derive the properties of the ionized gas even in the outer 
low surface brightness parts of the galaxies. However, the param- 
eters of the stellar populations are determined with much larger 
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uncertainties. We showed that analysis of azimutally averaged 
spectra at several de-projected galactocentric radii instead of the 
2D maps provides a more robust way to derive the radial de- 
pendencies of the stellar population properties in galaxies with 
well-defined axial symmetry. The main results of our study can 
be summarized as follows: 

- the six galaxies are quite massive with masses exceeding 2 x 
10 10 Mo; 

- the ionized gas and the stellar populations both indicate 
metallicities above the solar value; 

- five of the galaxies are currently forming stars at a rate of 
1-5 Mq, yr" 1 , which is typical for spiral galaxies at z - 0. 
The sixth galaxy shows no signs of star formation; 

- the five star- forming galaxies have mean mass-weighted stel- 
lar age ~ 5 Gyr and the passive one -12 Gyr; 

- four of the five star-forming galaxies show radial gradients 
of their ionized gas metallicity in the range from -0.022 to 
-0.058 dex kpc -1 . These values are typical for other spiral 
galaxies in the local universe. The fifth galaxy has a nearly 
uniformly distributed metallicity with a hint of a very low 
positive gradient of +0.007 dex kpc -1 ; 

- the radial dependence of the stellar population properties 
can be more robustly derived if azimutally averaged spec- 
tra at several de-projected galactocentric radii are analyzed. 
By this analysis we found indication of low negative radial 
metallicity gradients of the stellar populations in some galax- 
ies of up to -0.03 dex kpc -1 . Given the large uncertainties 
with which the stellar metallicity is estimated, the signifi- 
cance of these gradient is difficult to assess; 

- in the five star-forming galaxies the fraction of young stel- 
lar populations increases until 4-8 kpc and shows signs of a 
subsequent decrease. In four of them the fraction of old stars 
monotonically decreases in the disk and one galaxy shows a 
more complex behavior. 

- the passive galaxy has mostly old stars, with a possible small 
fraction of younger stars in the central few kpc; 

- the kinematic analysis indicates that the galaxies are relaxed 
systems that most likely have not experienced recent a major 
merger; 

- most of the SNe in our sample are projected on regions with 
metallicity and star formation rates above the galaxy aver- 
age, likely as a result of our target selection criteria and the 
radial metallicity gradients; 

- the BPT diagnostic diagram revealed that two of the galax- 
ies host AGNs. Another galaxy is on the border between the 
AGNs and the star-forming galaxies. Our analysis shows that 
the AGNs are not strong enough to affect the quantities de- 
rived from the total galaxy spectra. This AGNs may not be 
recognized in studies of host galaxies of SNe la at high red- 
shift; 

- the correlation of the SALT2 x\ parameter of the SNe and 
the Hubble residuals with the various parameters of the host 
galaxies did not lead to conclusive results. The low statis- 
tics and the small ranges spanned by the galaxy parameters 
render such an attempt still premature. We also note that the 
HRs of our SNe are on average positive, although their host 
galaxies have masses in the range where the other studies 
have shown negative HRs. 

In conclusion, we have demonstrated the viability to study 
the host galaxies of SNe la at low redshift using wide-field IFU 
spectroscopy at 4m-class telescopes. Intermediate-resolution 
spectra with sufficient S/N can be obtained out to the outer 
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low-surface brightness parts of the galaxies with a reasonably 
long exposure time of ~ 1.5 hours. Compared to a integrated 
spectroscopy and analysis of multi-color broad-band imaging, 
the IFU spectroscopy provides much more detailed information 
about the properties of the galaxies, e.g. metallicity and age gra- 
dients, detailed star formation histories, etc. In principle, the S/N 
of our data is sufficient to perform correlation analyses between 
the SN properties and the properties of the host galaxies at the 
location of the SN. However, our current sample it too small and 
suffers from strong selection biases to provide robust correla- 
tion results. The ongoing CALIFA survey may soon provide IFU 
spectroscopy of a larger sample of SNe la host galaxies, which 
will be a solid basis to further explore this path to study the 
SNe la progenitors and improve SNe la as distance indicators. 
We have tested the methodology and developed semi-automated 
tools that will allow expanding our analysis once the CALIFA 
data become available. 
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Appendix A: PMAS/PPAK instument details 

In this section we provide some details on the PPAK instrument 
that are relevant f or the data reduction. For a full descriptio n of 
the instrument see Verheijen et al. ( 2004| ), Roth et al. (2005 ), and 
|Kelz etal.| ( |2006l ). 

The PPAK fiber bundle consists of 382 fibers with 
2.7"diameter each, 331 of which (science fibers) are ordered into 
a single hexagonal bundle that covers a field-of-view (FOV) of 
72"x64". Thirty-six additional fibers form six mini-bundles of 
six fibers each (sky-bundles). The sky-bundles are evenly dis- 
tributed along a circle of rad ius ~ 90" and fa ce the edges of the 
central hexagon (see Fig. 5 in |Kelz et aL 2006| ). The remaining 15 
fibers are used for calibration and can only be illuminated with 
the PMAS internal calibration unit. 

The fibers are ordered into 12 slitlets. Each slitlet typically 
holds 28 science fibers, three sky-fibers from three different sky- 
bundles and one calibration fiber. The three sky-fibers are evenly 
distributed between the science fibers in the slitlet. The 12 slitlets 
are ordered to form a pseudo-slit. When projected onto the CCD 
in 2x2 binned mode, the separation between the spectra along 
the cross-dispersion direction is ~4.8 pixels and the full- width 
at half-maximum (FWHM) of the spectral traces is ~ 2.5 pixels. 
Between the slitlets there are gaps about two fibers wide, with 
the exception of the two central slitlets, which are separated by 
w ider gap of ~ 10 fi bers (see Fig. |B.l| and also Figs. 13 and 14 
in |Kelzetal.1 ( |2006l )). 

Because the fibers are circular, only -65% of the FOV is 
spectroscopically sampled in a single exposure. Therefore, at 
least three suitably offset pointings are needed to spectroscop- 
ically sample every point in the PPAK FOV. 



Appendix B: Data reduction 

The pre-reduction of the CCD images was performed with IRAF 
and the remaining reduction with our own programs written in 
IDL. 



B.1. Pre-reduction 

The PMAS 4kx4k CCD is read by four amplifiers and for each 
exposure four separate FITS files are created. These were in- 
dividually pre-reduced. The bias frames showed no significant 
large-scale structures and the bias was corrected by subtracting 
the average value computed from the CCD overscans. The im- 
ages were then trimmed, converted from ADUs into electrons 
using the gain values measured during the commissioning of the 
CCD and finally combined into a single 2D image, on which the 
spectra were oriented roughly along the rows. 

To create a master flat-field image, all halogen lamp expo- 
sures obtained during the run where summed. The separation 
between the spectra along the cross-dispersion direction is ~4.8 
pixels and the full-width at half-maximum (FWHM) of the spec- 
tral traces is ~ 2.5 pixels. As a result, the intensity of the pixels 
between the spectra is -20% of the peak value (e.g, see Fig. 18 in 
Kelz et al. 2006 ). Therefore, when all the images of the halogen 
lamp spectra are combined, there will be enough counts even in 
between the spectra. Moreover, because of the instrument flexure 
the positions of the spectra change by up to ~ 1-2 pixels, de- 
pending on the pointing of the telescope. As a result, in the com- 
bined image the intensities between the spectra were ~ 80% of 
the peak values. Because the intensity changes smoothly along 
the dispersion axis, each image row was smoothed with running 



median in a window of 20 pixels to normalize the master flat- 
field image. The flat-field image was divided into its smoothed 
version to derive the final normalized flat-field that contained 
only the pixel-to-pixel variations. 

Because of the relatively long exposures used, the images 
contain many pixels that are affected by cosmic ray (CR) hits. 
The complexity of PPAK images makes it very difficult to use 
conventional methods for CRs rejection. We experimented with 
different approaches and found that the CRREJ algorithm (de- 
veloped for WFPC on HS 70) Rave the best results. We tried dif- 
ferent parameter settings, each time carefully examining the CR 
masks to verify that no sky lines or galaxy emission lines were 
identified as CRs and at the same time as many as possible CRs 
were identified. It was occasionally necessary to manually mark 
some CRs. Finally, the values of the CR hit pixels were replaced 
with the values derived by interpolating the adjacent good pix- 
els. We used cubic-spline interpolation and the direction of the 
interpolation depended on the extent of the CR-affected region 
along the image rows and columns 

B.2. Spectra tracing 

The continuum lamp exposures that were obtained before and/or 
after the science exposures were used to trace the positions of the 
spectra on the CCD. The spectra tracing was performed in two 
steps. First, 20 columns at the middle of the image were aver- 
aged and used to determine the location of all 382 spectra. Then 
moving left and right and averaging 20 columns, the position 
of each spectrum at each CCD column was determined using the 
peak positions of the previous column as starting points. The po- 
sitions of the maxima were determined by fitting a parabola to ei- 
ther the three or the four highest intensity points. When the ratio 
between the two highest points was lower than 0.9, three points 
were used, in which case the parabola passes exactly through all 
three points and the maximum can be computed. Otherwise, a 
least-squares fit was used to determine the maximum. 

The space between the adjacent spectra is -4.8 pixels, which 
prevented us from using a more accurate scheme to compute the 
positions of the spectra. The traces determined at the first step 
were clearly oscillating with an amplitude -1/3 pixels and it was 
necessary to additionally smooth them. The traces were fit with 
a sixth-order polynomial function and the results were stored in 
a FITS file. 

Because of the mechanical flexures of PMAS the positions of 
the spectra moved slightly (at sub-pixel level) between the expo- 
sures, even within the sequence of exposures of a given galaxy. 
This effect can be accurately accounted for by a constant shift 
of the derived spectral traces. To derive these shifts, the average 
of 20 columns taken at the middle of the 2D image was cross- 
correlated with the same average in the halogen 2D image that 
was used to trace the spectra. In the subsequent reduction steps 
these shifts were always accounted for. 



B.3. Scattered light 

A small amount of scattered light is present in the images and 
it needs to be subtracted to achieve an accurate flux calibra- 
tion. To derive a 2D model of the scattered light, we used the 
gaps between the slitlets (Fig. B.l ), which contain no light from 
the object. Bands of five pixels width were extracted from the 



6 We use the IDL implementation available in the IDL Astronomy 
User's Library maintained by Wayne Landsman and available at http : 
//idlastro . gsf c . nasa . gov/ 
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Fig. B.l. Part of one of the science images that demonstrates the 
configuration of the fibers on the CCD. The calibration fibers 
are not visible because they are not illuminated during the sci- 
ence exposures. The calibration fibers are normally last within 
the slitlets. Slitlet #6 holds two calibration fibers, and two more 
are placed at the beginning of the fiber sequence and one is at 
the end. 



gaps (wider band was extracted between the sixth and the sev- 
enth slitlets) using the information from the spectra tracing. The 
bands were converted into ID vectors by running median in 5x5 
pixel window and the positions along the y-axis were also com- 
puted. The scattered light pattern along the dispersion axis is 
fairly complex and the attempts to directly fit a 2D polynomial 
surface yielded no satisfactory results. Instead, the scattered light 
vectors were cubic- spline-interpolated along all image columns 
(in the spatial direction) to derive the amount of scattered light 
at each image pixel. The resulting 2D model image of the scat- 
tered light was smoothed by 50x50 pixel running median and 
subtracted from the original image. 

B.4. Spectra extraction and cross-talk correction 

The spectra were extracted from the 2D images using a simpli- 
fied version of the optimal extraction algorithm of Home ( |1986| ). 
In this method the spatial profile of the spectral trace is deter- 
mined from the data themselves and is used to form a weighted 
sum of the pixels in a given aperture, giving larger weights to the 
pixels that received more light. In the conventional longslit spec- 
troscopy the profile is determined by the seeing and the telescope 
guiding, and may change considerably from exposure to expo- 
sure. In the case of PPAK the spatial profiles are almost entirely 
determined by the size of the fibers and have a stable shape that 
can be accurately described by a Gaussian. Thus, we determined 
the FWHMs of the Gaussians that describe the profiles of all 
382 spectra and used this information to compute the weights for 
the optimal extraction. Series of halogen lamp spectra were ob- 
tained and added to increase the signal-to-noise ratio (S/N). For 
each image column 382 Gaussians were simultaneously fitted to 
derive the FWHMs. The fitting was performed with the MPFIT 



IDL packag^] and special care was taken to provide accurate 
starting values; otherwise the fit with so many Gaussians would 
fail. The resulting 2D FWHM image (382 spectra by 2048 spec- 
tral elements) was somewhat noisy and was smoothed by run- 
ning median within 11x11 elements window. The mean FWHM 
is 2.46 pixels and the covered range is 2.12-2.85. 

With the positions of the spectra and the spatial profiles de- 
termined, the spectra were extracted with the optimal extraction 
algorithm. The separation between the spectra within the slitlets 
is fairly constant ~ 4.8 - 4.9 pixels, both from spectrum to spec- 
trum and along the trace. The extraction apertures were auto- 
matically computed as half the average distance from closest 
spectrum. This means that for a given spectrum the extraction 
aperture may be slightly asymmetric with respect to the trace 
center and that different spectra may have slightly different aper- 
ture sizes. 

Within a given slitlet the spectra are separated by only ~ 
2xFWHM, which leads to an effect known a s cros s-talk, i.e. 
contamination by adjacent spectra^ Sanchez (2006) proposed 



the iterative Gaussian-suppression technique and we followed 
this approach with small modifications. After the first optimal 
extraction, the contamination from the adjacent spectra for a 
given spectrum was estimated as follows. Using the estimated 
flux, trace positions and FWHMs of the Gaussians that describe 
the spatial profiles of the adjacent spectra, and taking into ac- 
count the finite size of the extraction apertures, we computed 
the height of the Gaussian functions that would have produced 
the fluxes of the adjacent spectra were. This was performed for 
all columns and the contamination of the adjacent spectra was 
subtracted. Then a new optimal extraction was performed on the 
corrected image using the same extraction aperture as before. 
This was performed for all spectra and the process was iterated 
once more. 



B.5. Wavelength calibration 

The HgNe and ThAr arc-lamp spectra were used to derive sepa- 
rate wavelength calibrations for each individual spectrum. The 
arc-lamp spectra were extracted in exactly the same way as 
the science ones and the positions of all Hg, Ne and CcQ lines 
with known laboratory wavelengths were determined by a least- 
squares fit with Gaussian function. A sixth-order polynomial 
fit was then used to derive accurate pixel-to-wavelength trans- 
formations. The spectra extend from about 3700 A to 7000 A. 
The bluest line in the HgNe spectrum is at ~4047 A. Thus, the 
wavelength solution below -4000 A may not be accurate. On the 
other hard, the two reddest arc lines are at 6717 A and 6929 A. 
Because of the field vignetting (see below) the 6929 A line was 
not detected in all fibers and as a result for those fibers the red 
part of the wavelength solution may also be inaccurate. To work 
around this problem, the 15 calibration fibers that were illumi- 
nated with a ThAr lamp were used to improve the wavelength so- 
lutions in the two extremes of the covered wavelength range. The 
dispersion over the wavelength range is between 1.5 A pixel -1 
and 1.7 A pixel -1 



7 available at |http : //cow .physics .wise . edu/~craigm/idl/ 
|idl.htmlj 

8 In the case of PPAK the effect is small but we nevertheless chose to 
correct for it. 

9 few Cadmium lines have also been identified in the spectra 
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B.6. Background subtraction 

The spectra extracted from fiber-fed or IFU spectrographs con- 
tain not only the object light, but also the sky background, which 
needs to be subtracted. For faint objects the background may ex- 
ceed the object flux by far and accurate background subtraction 
is essential to obtain reliable results. The PPAK design with the 
separate sky-fibers, which are distributed evenly between the sci- 
ence fibers, gives the flexibility to apply different approaches to 
the background subtraction. The following method was found to 
work best. The sky-fibers were first examined to check if some of 
them were accidentally contaminated by background objects. No 
such cases were found. Each element of the 331 science spectra 
was labeled with a pair of numbers (A t , N), its wavelength and 
fiber number as is recorded on the CCD, respectively. For a given 
element (Ai, N) all sky- spectra were interpolated at a common 
wavelength Ai using forth-order B -spline interpolation/^) which 
has one of the best interpolation properties of all known inter- 
polation schemes (e.g., |Thevenaz et al.|20 00). The interpolated 
sky-spectra were fitted as a function of their fiber position with 
least-squares fifth order polynomial fit and the value of the sky 
at the element position N was recorded. Repeating this for all 
elements gives a 2D model of the sky-background, which was 
subtracted from the science fibers. The sky- subtracted 2D im- 
ages were carefully examined and no visible systematic resid- 
uals were noticed, except at the position of the few strongest 
sky-lines. 



B.7. Flux calibration 

The flux calibration of the spectra was a two-step process: (i) 
normalization for the fiber-to-fiber throughput variations and 
(ii) correcting for the wavelength-dependent system response. 
To account for the fiber-to-fiber throughput variations the se- 
ries of twilight sky exposures were used. The individual images 
were added to increase the S/N and all spectra were extracted 
and wavelength-calibrated in the same way as the science ones. 
Unfortunately, with the new larger CCE^the images suffer from 
vignetting at the corners. The sensitivity of the parts of the spec- 
tra that are projected at the CCD corners f alls sh arply as function 
of the distance from the CCD center (Fig. |B.2| ). The relation be- 
tween the on-sky fiber positions and the fiber position on the 
CCD is such that the most affected fibers are t hose adjacent to 
the central hexagon of 37 fibers (see Fig. 7 in |Kelz et al.|2006| ) 
and not the outermost fibers. The fibers of the central hexagon 
are projected at the middle of the CCD and are practically un- 
affected by the vignetting. 

To derive the relative throughput of the fibers, all spectra 
were interpolated on a common wavelength frame. The 37 spec- 
tra of the central hexagon were averaged and all spectra were 
divided by this average spectrum. The resulting ratios were least- 
squares fitted with polynomials of different degrees (high orders 
were needed for the fibers most affected by the vignetting) to 
derive the throughput map of the fibers. The map is shown in 
Fig. B.2 and the vignetting at the corners is clearly visible. This 
map was used to equalize the throughput of the fibers by dividing 
all science exposures by it. 

During the analysis it was found that for several of the fibers 
that were most severely affected by the vignetting the drop of 
the throughput in the extreme blue and red parts could not be 
corrected well. This resulted in several spaxel spectra having in- 




10 forth-order B-spline was used for all spectra interpolations 

11 the detector was upgraded to 4kx4k CCD in October 2009 



Fig. B.2. Throughput map of the PPAK fibers. Note the vi- 
gnetting in the corners. 



correct relative flux calibration at the two extreme ends. This 
affected the spectrum fitting of these spaxels with STARLIGHT. 
Because the Ha flux was also affected, the estimation of the ex- 
tinction was unreliable and hence part of the emission line anal- 
ysis. The 'bad' spaxels were masked in the 2D maps of the af- 
fected properties. It should be noted, however, that our primary 
gas metallicity indicator, the 03N2 method, was not affected be- 
cause it only involves flux ratios of emission lines at close wave- 
lengths. 

The two spectrophotometric standard stars that were ob- 
served on 15 November were used to derive the sys- 
tem response. The spectra of the standards were ex- 
tracted, wavelength-calibrated, background- subtracted and fiber- 
throughput-corrected in the same way as the rest of the science 
exposures. The presence of gaps between the fibers in combi- 
nation with the point source nature of the standards makes it 
difficult to recover the whole flux of the star from a single expo- 
sure. In addition, because of the differential atmospheric refrac- 
tion (DAR) there may be significant differential flux losses. By 
taking the observations of the standards at five different points 
on the PPAK FOV, we aimed to minimize these effects. First, the 
differential flux losses may cancel out when the five spectra are 
summed. Second, it gives more control, for example, too discard 
a clearly deviating observation. 

To minimize the flux losses instead of summing all spec- 
tra that contain object flux, we adopted a different approach. 
The spectra were interpolated at a common wavelength scale 
and each fiber was associated with its on-sky position. For each 
wavelength plane a symmetric 2D Gaussian function (without 
constant term) was fitted and the flux was taken to be the in- 
tegral over the Gaussian. So extracted spectra were compared 
with the result of simply suming the fibers. The spectra extracted 
with the Gaussian fitting were clearly more consistent between 
the five pointings. With the exception of one pointing for each 
standard, the systematic difference between the spectra were less 
than ~ 3 - 4%. The deviating spectra for each star were dis- 
carded and the remaining four were averaged. After correcting 
for the atmospheric extinction using a mean extinction curve for 
Calar Alto Observatory, the ratio between the observed counts 
and the tabulated fluxes of the standard were fitted with a poly- 
nomial function to derive the sensitivity function of the system. 
The sensitivities derived from the two standards agreed well and 
the average of the two sensitivities was used to calibrate all sci- 
ence observations. 

B.8. Final data-cubes 

The final step in the reduction was to combine the three point- 
ings into a data cube. The spectra were re-sampled into a com- 
mon linear wavelength scale with a sampling of 1.5 A per el- 
ement. Taking into account the offset between the exposures, 
the circular fibers were mapped onto an output grid of square 
2 // x2 // pixels. The flux of a given fiber was distributed between 
the output pixels according to the overlap area between the fiber 
and the output pixels. No attempt was made to first shrink the 
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Fig.B.3. Example of the effect of the stellar 
spectrum subtraction in the regions of Yip and 
Ha emission lines. In the top panels we show 
the observed spectrum with the best fit over- 
plotted. The bottom panels show the observed 
spectrum minus the best fit. Clearly, the fluxes 
of the emission line cannot be accurately mea- 



4700 4800 4900 5000 51 00 6520 6560 6600 6640 sured without a proper subtraction of the stellar 

Rest wavelegth [A] Rest wavelegth [A] spectrum 



fibers and emulate the D rizzle algorithm used for HS T images 
(Fruchter & Hook 2002 ); with only three barely overlapping 
pointing this would not gain any improvement. However, with 
five or more carefully selected pointings it may be possible to 
use Drizzle to improve the spatial sampling of the final cubes. 

Two important effects were taken into account two when cre- 
ating the data-cubes - the non-photometric condition and the dif- 
ferential atmospheric refraction (DAR). To account for the first, 
for each galaxy and pointing the total flux within an aperture of 
12" centered on the galaxy core was computed and the observa- 
tions of given galaxy were scaled to the exposure with the high- 
est flux. With one exception all corrections were smaller than 
15%; the remaining pointing required a correction factor of ~2. 

The observations were obtained at an airmass of less that 1.3 
and the effect of DAR was small, but we nevertheless corrected 
for it. To compute the DAR we used the formulas from Szokoly 
2005) and the ambient conditions as recorded in the FITS image 



headers. To combine the three pointings, we shifted the positions 
of the fibers for each wavelength slice along the x and y axes by 
the values computed from the value of the DAR at that wave- 
length (with respect to the refraction at 6300 A) and the effec- 



tive parallactic angle of the observations (see also Walsh & Roy 



11990] ). For both the DAR-corrected and un-corrected data-cubes, 
the positions of the galaxy nuclei were computed in the blue 
(average in 3800-4000AA) and red (6600-6800AA) wavelength 
ranges. The two positions measured in the DAR-corrected cubes 
matched very well, but not in the un-corrected cubes, which 
shows that the DAR was correctly accounted for. 

Three of the galaxies in our sample also have SDSS spec- 
troscopy. This allowed us to check the relative flux calibration 
of our spectra. To emulate the SDSS spectroscopy, we extracted 
from our data-cubes the flux within an aperture of 3" diameter 
centered on the galaxy nucleus and scaled it to match the mean 
flux level of the SDSS spectra. The comparison between the 
SDSS spectra and those extracted from our date-cubes is shown 
in Fig.[T] It shows that the relative flux calibration of our spectra 
is excellent and matches SDSS to within a few percent. 

To set the absolute flux scale of the data-cubes we used the 
SDSS imaging. SDSS g and r magnitudes of the galaxies were 
computed within an aperture of 20" diameter. Spectra within the 
same aperture size were extracted from the data-cubes and syn- 
thetic g and r magnitudes were computed. The g and r scale 
factors that made the synthetic magnitudes match the observed 
ones were computed and the average of the two was applied to 
the data-cubes. It should be noted that the g and r scale factors 
coincided to within a few percent, which supports our conclusion 
that the relative flux calibration is accurate. 



Appendix C: Measurement of emission line fluxes 

Five of the six galaxies show strong nebular emission lines. 
This enabled us to study the properties of the gas phase in the 
galaxies through measuring the fluxes of the most prominent 
emission lines [On]/l3727, H/5, [O m] /U4959/5007, Ha, and 
[Nil] /U6549/6584, and when possible also [Sn] /U67 16/6731. 
In spectra of galaxies the nebular emission lines are super- 
imposed on the underlying stellar continuum. The stellar ab- 
sorption lines can bias the measurement of the emission line 
fluxes, an effect that is especially prominent in UJ3 (Fig. |B.3 ). 
Therefore, the stellar continuum needs to be subtracted first to 
measure the emission line fluxes accurately. For this we used 
the STARLIGHT software ( |Cid Fernandes et al.|2005) . Briefly, 
after masking the regions of known nebular emission lines, tel- 
luric absorptions, and strong night-sky emission lines, the spec- 
trum was fitted with a linear combination of model spectra of 
single stellar populations (SSP) of different ages and metallici- 
ties. The best fit was subtracted to derive the pure emission line 
spectrum, where the emission line fluxes can be measured with- 
out bias, see Fig. B.3 The parameters of the fitted SSPs can be 
used to derive the properties of the galaxy stellar population, 
but here we regarded the fits as a means to subtract the under- 
lying stellar absorption. All spectra that had S/N higher than 5 
at ~ 4600A were fitted and the emission line fluxes were mea- 
sured on the continuum-subtracted spectrum. For the remaining 
spectra the measurements were made without continuum sub- 
traction. A weighted non-linear least-squares fit with a single 
Gaussian plus a linear term was performed for each emission 
line, and the area below the Gaussian was taken as an estimate 
of the flux. The weights were the error spectrum produced during 
the spectral extraction. Including the linear term is necessary to 
account for possible systematic residuals due to imperfect sub- 
traction of the the stellar absorption spectrum. The close lines 
Ha, [Nn]/i6549 and [Nn]/i6584 were simultaneously fitted. 
[On] ,13727 is a blend of two lines - [On] ^3726.04/3728.80. 
This blend was fitted with two Gaussians with equal width and 
a fixed separation of 2.76 A. In addition, a robust standard de- 
viation of the adjacent continuum was also measured from the 
regions more then 3<x away from the line centers. 

The uncertainty of the line fluxes was estimated by prop- 
agating the uncertainties of the fitted amplitude and cr of the 
Gaussians. To check the reliability of this estimate, we per- 
formed Monte Carlo simulations. We assumed that the lines have 
a Gaussian shape with cr = 2, representative for our spectra, 
and that they were superimposed on constant background with 
standard deviation one and mean zero. The noise across the line 
was assumed to follow a Poisson distribution - if a given da- 
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turn has N counts, the l<x uncertainty of this datum is (1 + N) l/2 . 
This noise model was used to generate pseudo-random numbers, 
which were added to the line and the line was fitted (weighted 
by the noise model). This was repeated 1000 times and the mean 
values and the standard deviation of <x, the amplitude, and the 
area below the Gaussian were computed. This procedure was re- 
peated for Gaussians with amplitudes 2 to 1000, i.e. exceeding 
the continuum noise by the same numbers. 

The results of these simulations show that propagating the 
errors of the fitted parameters to the total flux overestimates the 
uncertainty by about a factor of 2. The simulations also showed 
that the area below the Gaussian can be recovered with ~ 30% 
accuracy and without significant bias for lines that exceed the 
continuum noise by only three times. The S/N of the measured 
line flux was tabulated as a function of the ratio of the fitted 
amplitude of the Gaussian to the standard deviation of the adja- 
cent continuum and was used to assign realistic errors to the line 
fluxes measured by fitting a Gaussian function. For example, a 
10% accuracy is achieved when the amplitude of the Gaussian 
is 20 times the continuum noise. A possible concern is that the 
shape of the lines in the real spectra may deviate from Gaussian, 
in which case additional uncertainty will be introduced. Close in- 
spection of the line fits revealed that the line shape in our spectra 
is well represented by a Gaussian function and we do not expect 
problems related to a non-Gaussian line shape. 



Appendix D: Online figures 
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Fig. 7. From left to right: SDSS color images of the galaxies, the light- weighted average stellar population age (log and metal- 
licity (Z*)l, and the visual extinction Ay determined by STARLIGHT fits to the stellar spectra. 
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Fig. 8. STARLIGHT fits to the total galaxy spectra. For each galaxy we also show the ratio between the observed and the best-fit 
spectrum. The vertical brown lines in the spectrum panels show the location of the strongest night-sky lines, which could not be 
cleanly subtracted and the corresponding wavelength regions were excluded from the fits. To the right of the plots are shown the 
population vectors and the mass-fraction vectors along with the extinction, the mass- and light- averaged age and metallicity, and the 
contribution of young (age < 300 Myr), intermediate (300 Myr < age < 2.4 Gyr) and old (age > 2.4 Gyr) stellar populations. The 
short brown bars show the ages of the SSPs used in the fits and the two vertical dotted lines separate the young, intermediate and 
old populations. 
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Fig. 9. Mass- and light-weighted average stellar populations metallicity (left panel) and age (right panel) from the STARLIGHT fits 
as function of the de-projected galactocentric distance. The small black dots are the measurements on the individual spaxels and the 
large blue dots are from the azimuthally averaged spectra. 
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Fig. 12. Same as in Fig.[i] but for the azimutally averaged spectra. 
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Fig. 13. From left to right: SDSS color images of the galaxies and the maps of the fraction of young (age < 300 Myr), intermediate 
(300 Myr < age < 2.4 Gyr) and old (age < 2.4 Gyr) stellar populations. 
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Fig. 14. Compressed population vectors corre- 
sponding to the contribution of young (age < 
300 Myr), intermediate (300 Myr < age < 2.4 
Gyr), and old (age > 2.4 Gyr) stellar popula- 
tion to the formation of the observed spectra as 
function of the de-projected galactocentric dis- 
tance. The small black dots are the measure- 
ments obtained from the individual spaxels and 
the large blue dots are from the azimuthally av- 
eraged spectra. 
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Fig. 15. From left to right for each galaxy we show the color SDSS image, Ha velocity map, the star velocity map, and the difference 
between them. The x, y coordinates are in arcsec with respect to the map centers. The orientation of the images is north - up, east - 
left. 



